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Abstract 

We present complete results of a next-to-leading order calculation of the 
production of gaugino-like charginos (x^) and neutralinos (x^) in association 
with gluinos (g) at hadron colliders, including the strong corrections from the 
exchange of colored particles and sparticles. Adopting a variety of models for 
the sparticle mass spectrum, including typical supergravity (SUGRA) models 
and a light gluino model, we provide predictions for total and differential cross 
sections at the energies of the Fermilab Tevatron and CERN Large Hadron 
Collider (LHC). 
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I. INTRODUCTION 



Weak scale supersymmetry (SUSY) |l| is a theoretically attractive extension of the Stan- 
dard Model of particle physics. Supersymmetric theories can solve the Higgs hierarchy puzzle 
0, break electroweak symmetry radiatively at low energies [Q, and explain the unification 
of the gauge couplings at a high energy scale . SUSY introduces a superpartner for each 
Standard Model particle with the same quantum numbers, but for a difference in spin of 
1/2. If supersymmetry were exact, these superparticles would be degenerate in mass with 
their Standard Model partners. However, SUSY can be broken softly in such a way that its 
attractive features survive, while the superpartners become heavy enough to evade current 
limits from collider searches |^ . A powerful and general parametrization for the soft SUSY- 
breaking terms is provided by the Minimal Supersymmetric Standard Model (MSSM) 0. 
Soft breaking mechanisms require that the superpartners remain lighter than a few TeV, 
and thus there is reason to expect that high energy investigations, such as those at LEP II, 
Run II of the Fermilab Tevatron, and the CERN Large Hadron Collider (LHC) will discover 
supersymmetric particles or place strong constraints on supersymmetric models. 

The search for direct experimental evidence of supersymmetry at colliders requires a 
good understanding of theoretical predictions of the cross sections for production of the su- 
perparticles. In the case of hadron colliders, where collisions of strongly interacting hadrons 
are studied, the large strong coupling strength (as) leads to potentially large contributions 
beyond the leading order (LO) in a perturbation series expansion of the cross section. To 
have accurate theoretical estimates of production rates and differential cross sections, it is 
necessary to include corrections at next-to-leading order (NLO) or beyond. Next-to-leading 
order calculations of the hadroproduction of gluinos and squarks , top squarks [H , sleptons 
PJlOd, and gauginos [|10[ have been published, including our brief report on associated pro- 



duction of gauginos and gluinos [|Tl|]. In this paper, we provide a detailed exposition of our 
calculation of associated production, and we present new predictions of total and differential 
cross sections for a variety of assumptions about the superpartner mass spectrum. 
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Associated production of a gaugino (x) with a gluino (g) or with a squark (g) is po- 
tentially a very important production mechanism. Associated production processes are 
semi-weak in that they involve one somewhat smaller coupling constant than the produc- 
tion of a pair of colored sparticles. However, in popular models of SUSY breaking |T2|,|13[, 
the mass spectrum favors much lighter masses for the color-neutral, low-lying neutralinos 
and charginos than for the colored squarks and gluinos. Their lighter mass means that the 
phase space for production of neutralinos and charginos, and the relevant partonic lumi- 
nosities, will be greater than that for gluinos and squarks. This effect is potentially decisive 
at a collider with limited energy, such as the Tevatron. Indeed, as our numerical results 
show, the extra phase space may more than offset the smaller coupling, and gauginos may 
be produced more copiously than squarks at the Tevatron. Another point in favor of asso- 
ciated production is the relative simplicity of the final state. For example, the lowest lying 
neutralino is the (stable) lightest supersymmetric particle (LSP) in supergravity (SUGRA) 



models 12|, manifest only as missing energy in the events, and it is the second lightest in 



gauge-mediated models [|T3|. The charginos and higher mass neutralinos may decay leptoni- 



cally leaving a lepton signature plus missing transverse energy; relatively clean events ensue. 
Furthermore, associated production may be the best channel for measurement of the gluino 



mass 
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In this paper we present our complete NLO calculation in SUSY-QCD of the hadropro- 
duction of a gluino in association with a gaugino, including contributions from virtual loops 
of colored sparticles and particles, and three-particle final states in which light quarks or 
gluons are emitted. We extract the ultraviolet, infrared, and coUinear divergences by use 
of dimensional regularization and illustrate how they may be absorbed by the usual renor- 
malization and mass factorization procedures. In computing the virtual contributions, we 
encountered divergent four-point functions that had not been evaluated previously. We use 
a combined analytic and numerical phase space slicing method to treat the contributions 
from real emission of light particles. Associated production was calculated at LO some years 



ago |114| , p!5| . Our reason to focus first on the x plus g final state, rather than the x plus q 



final state, is that the LO cross sections for x + g are 3 to 6 times greater than those for 
X + at the energy of the Tevatron when the mass rrig = rriq = 300 GeV, and 6 to 15 times 
greater when rrig = niq = 600 GeV. These comparisons are pertinent for the hghter mass 
neutralinos Xi2 chargino x^- obtaining the q cross sections, we sum over five flavors 
of degenerate squarks and antisquarks. 

Our analysis is general in that it is not tied to a particular SUSY breaking model. We 
can provide cross sections for arbitrary gluino and gaugino masses, and, indeed, the values 
of the cross sections at Tevatron and LHC energies depend crucially on the sparticle masses. 
Mass generation in supersymmetry is accomplished in a hidden sector and transmitted to 
the MSSM fields. In SUGRA models [|12[, transmission is through gravitational interactions. 



while in gauge-mediated [|T3[ and gaugino-mediated |ll6|-[T^ models, it occurs through gauge 



interactions. Anomaly mediated SUSY breaking, also gravitational in origin, is a fourth 
possibility [|19|. In the presentation of predictions for cross sections, we consider illustrative 



mass spectra typical of each scenario and consistent with bounds established from current 
data p(]| , pT|] . We also examine the phenomenologically open case of a gluino with mass 
light compared to the SUSY scale, rrig ~ 30 GeV. This possibility arises in some models of 
gauge-mediated SUSY breaking E^l 



As is shown in detail below, the NLO corrections to associated production are generally 
positive, but they can be modest in size, ranging in the SUGRA model from a few percent 
at the energy of the Fermilab Tevatron to 100% at the energy of the LHC, depending 
on the sparticle masses. In the light-gluino case, NLO contributions increase the cross 
section by factors of 1.3 to 1.4 at the energy of the Tevatron and by factors of 2 to 3.5 at 
the energy of the LHC. Owing to these enhancements, collider searches for signatures of 
associated production will generally discover or exclude sparticles with masses larger than 
one would estimate based on LO production rates alone. More significant from the viewpoint 
of reliability, the renormalization and factorization scale dependence of the cross sections is 
reduced by a factor of more than two when NLO contributions are included. 

At Run II of the Fermilab Tevatron, for an integrated luminosity of 2 fb~^, we expect 



that 10 or more events could be produced in each of the hghter gaugino channels of the 
SUGRA model, gxi, 9X2^ a-iid gxf, provided that the gluino mass rUg is less than 450 GeV. 
The cross sections for the three heavier gaugino channels, ^Xs; 9x1^ ^ind gx^, are smaller 
by an order of magnitude or more than those of the lighter gaugino channels. In the light 
gluino model that we consider, more than 100 events could be produced in the three lighter 
gaugino channels provided that the common GUT-scale fermion mass mi/2 is less than 400 
GeV, and as many as 10 events in the three heavier gaugino channels as long as mi/2 is 
less than 200 GeV. At the higher energy and luminosity of the LHC, at least a few events 
should be produced in every channel in the SUGRA model and many more in the light 
gluino model. 

The shapes of the rapidity distributions of the gauginos are not altered appreciably by 
NLO contributions, but the locations of the maximum cross section in transverse momen- 
tum (pt) are shifted to smaller values by NLO contributions. At LHC energies where the 
contribution of the qg initial state is important, modifications of the px spectra can be 
pronounced. 

We begin in Sec. y with a brief review of the LO calculation in order to introduce 
our notation. In this section, we also introduce the SUSY breaking models we adopt and 
summarize salient aspects of their predicted mass spectra. This discussion is followed in 
Sec. m by a detailed presentation of our NLO (9(a;a;|) calculation. We present partonic 
scaling functions in Sec. |^ as well as predictions for inclusive and differential cross sections 
at Tevatron and LHC energies. Our conclusions may be found in Sec. 0. In Appendices - 
0, we present detailed analytic results related to the NLO calculation. 

II. LEADING ORDER PRODUCTION OF GAUGINOS AND GLUINOS 

We begin with the Born level cross sections for the partonic processes 

^ gxtQQ ^ 9xt^ (i) 
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derived first in |TJ,|T5|. In anticipation of tlie renormalization and mass factorization of 
the NLO contributions, we proceed in the n = 4 — 2e dimensions of standard dimensional 
regularization. We assume that there is no mixing between squarks of different generations 
and that the squark mass eigenstates are ahgned with the squark chirality states, equivalent 
to the assumption that the two squarks of a given flavor are degenerate in mass. We ignore 
the n/ = 5 light quark masses in all of the kinematics and couplings, and thus study the 
production of gaugino-like charginos and neutralinos, but not the production of Higgsino-like 
ones ||2^. We assume further that the entries in the chargino and neutralino mass matrices 
are real, and thus that the unitary transformations from the {B, W^, Hi, H2) and {W^ , 
H^) bases to the (x?, X2; X4) and (x^, xt) bases are given by orthogonal matrices. A 
result of this convention is that it is possible for the mass of one or more neutralinos to be 
negative inside a polarization sum. The chargino masses are chosen to be positive as may 
be done for Dirac fermions. 

The Dirac matrix 75 (or equivalently the projectors Pl{r) = (1 ='=75)/2) appearing in the 
gaugino and gluino couplings is treated in the 'naive' scheme in which it anti-commutes with 
all of the other 7^ matrices. This scheme is acceptable at the one loop level for calculations 
free from anomaly. In evaluating the Feynman diagrams involving Majorana and explicitly 
charge-conjugated fermions, we have followed the approach described in Ref. [p5| . 

We express our leading order results in terms of the Mandelstam variables 

S = {Pa+Pl,f (2) 

t = (Pb- Pif h = {Pb - P2Y - fni"^ t2 = {Pb - P2f - 

U = {Pa-P2f Ui = {pa - P2f - ^2 = (Pa " ^2)^ - ^T^2^ 

where Pa-, Pb, Pi, and p2 refer to the four-momenta of the incoming quark, the incoming anti- 
quark, the produced gluino, and the produced gaugino, respectively. Variable mi denotes 
the mass of the gluino and m2 that of the gaugino. The incoming partons are treated as 
massless. The momenta are on mass-shell, Pa = Pb — O5 Pi — ^ly and P2 = m\. The 
invariants obey the relation s + t + u = m\ + m^. 
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After the n-dimensional phase space integration we obtain the lowest order partonic 
differential cross section, 



dt2 du2 r(l — e) 



t2 U2 — mis 

fi"^ s 



Q{t2U2- ni\s) (3) 



xe(s - (mi + ) 5{s + t + u - m{ - m^) \M 



where = (47r) The arbitrary scale fi is introduced, as usual, to provide the correct 

2 

mass dimension for the coupling in n dimensions; is the leading order matrix element 

summed over the colors and helicities of all of the outgoing particles, and averaged over the 
colors and helicities of the incoming ones. The indices label the incident partons. For 
neutralino production at leading order, the partons are quarks and antiquarks of the same 
flavor. For chargino production, the incident quarks and antiquarks have different flavor. 
At next-to-leading order, quark gluon initial states contribute also. 

As is shown in Fig. |l], the Born matrix element for associated production of gluinos 
and gauginos proceeds via t— or u— channel exchange of a squark. In the case of charged 
gauginos, only the left-handed chiral squarks participate, whereas neutral gauginos receive 
contributions from both left- and right-handed chiral squarks. Furthermore, in the case of 
charged gaugino production, the squarks exchanged in the t— and u— channels correspond 
to different flavors, while in the neutral case the t— and u— channel squarks have the 
same flavor. Under our assumption that the squark mass eigenstates correspond to squarks 
of definite chirality, the (massless) incoming quark and anti-quark are forced to have a 
particular helicity, and thus the sets of graphs in which a right-handed chiral squark is 
exchanged cannot interfere with those mediated by a left-handed chiral squark. The matrix 
element has the analytic form [|l^ 
Svr as 



\M 



B\ 



9 



Xttit2 2 XtuS 7711 7712 ^ X„ Mi M2 



(4) 



(t-m|)2 (t_m2j(M-m?J (M-m?J2 

where as = gl/'^T^ is the coupling between gluinos, squarks, and quarks; m^^^^j is the mass 
of the squark exchanged in the t — (m— ) channel graph; and the X represent the gaugino 
interactions with quark and squark. 



For production of a neutralino of type x", the X are [p6 



Xt = X, 



tu 



Xu — 2 



(5) 



sin 6*14^ cos6'vi/ 

In the expressions above, e is the electron charge, 6w the weak mixing angle, Tg the third 
component of the weak isospin for the squark, and Cg is the charge of the quark in units of 
e. For up-type quarks Cq = 2/3 and for down-type quarks Cq = —1/3. The matrix A^' is the 
transformation from the interaction to mass eigenbasis defined in Ref. ||26|. The expressions 
for production of a positive chargino of type xt 



Xt 



sin^ e 



\Vk^ 



tu 



W 



sin^ e 



Re (y^i X, 



w 



. 2a l^'^il 
sm Uw 



(6) 



For the negative chargino Xk they have the form. 



Xt 



sin^ e 



\Uki\ 



X, 



tu 



W 



sin^ e 



ReiV,\Uki), Xu 



w 



sm 6w 



(7) 



Matrices U and V are the chargino transformation matrices from interaction to mass eigen- 



states defined in Ref. 26 



To compute cross sections for hadropro duct ion. 



hah^ g Xk X, hah^ g Xk X, (8) 

where ha and hi, label the incoming hadrons, one must convolve the partonic cross section 
with the parton distribution functions. In the high energy scattering limit, one may neglect 
the mass of the incoming hadrons compared with their momenta, and obtain 



S={Pa + Pbf 

T={Pt- P2f Ti = (P5 - P2f - T2 = {Pb - P2f - 

U={Pa- P2f U, = {Pa - P2f ' U2 = {Pa " - 



(9) 



in which Pi indicates the momentum of hadron h^. We define Xi by the relations 



S = XaXbS, t2 = Xi,T2, U2 = XaU2. 



(10) 
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The convolution with the parton distribution functions may be written ||^ as 

(S,T,[/,/i^)= (11) 



dT2 dU2 



^ [ dXa I dXb Xaf^" {Xa, Xbfj'' (Xfc, ^ ^'J' . 

In this equation fip refers to the factorization scale, and d'^&ij/ dt2du2 is the hard cross 
section, equal to the Born cross section at leading order. The lower limits of integration on 
the convolution are 

= S + fTa' ^^^^ 
_ -Xg U2-ml + m\ 

The differential cross section in the transverse momentum (p-r) and rapidity (y) of the 
gaugino is related to the differential cross section in U2 and T2 by 

rfpr dy dT2 dU2 ' 

with 

2 T2U2 2 ^2^2 2 

Pt = ^2 = "^2> (15) 

b s 

. = ilog(|). (16) 

The total cross section is obtained by integrating over the full range of transverse momentum 
and rapidity, 

a{S,fip)= dpT dy — (17) 

Jo J -y^^^{pT) apTay 

d^^y ■ 



The limits of integration are 



and 



pT''(y) = — ^ J (S + ml- m\f -AmlS Cosh2 (?/), (19) 

^ ^' 2v^Cosh(y)V^ 2 II 2 \yi^ K ) 



y'^'^^PT) = ArcCosh I ^ ^ : 1 . (20) 

' 2 {p\ + mi) ' ^ ^ 
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A. Super symmetry Breaking Models 



The physical gluino and gaugino masses that we use, as well as the gaugino mixing 
matrices, are based on four popular SUSY breaking models plus a fifth scenario in which 
the gluino mass is relatively light. 

For our default minimal SUGRA scenario |T2], we select the common scalar and fermion 
masses at the GUT scale to be mo = 100 GeV and mi/2 = 150 GeV. The trilinear coupling 
Aq = 300 GeV, and the ratio of the Higgs vacuum expectation values, tan/5 = 4. The 
absolute value of the Higgs mass parameter /z is fixed by electroweak symmetry breaking, 
and we choose /i > 0. (Our sign convention for Aq is opposite to that in the ISASUGRA 
code |2^.) For this scenario, the neutralino masses m^o ^ are 55, 104, 283, and 309 GeV with 
m^o < inside a polarization sum. The chargino masses i^^f^ are 101 and 308 GeV and 
therefore almost degenerate with the masses of X2 4- The gluino mass rrig is 410 GeV, and 
the squark mass is 359 GeV. All of these masses are above the exclusion limits established 
from LEP and Tevatron collider data [^,^. Since the gluino and gaugino masses vary 
principally with mi/2, we freeze the values of the other four parameters, and we vary mi/2 
over the range 100 to 400 GeV. The squark, gluino, and gaugino masses all increase as mi/2 
increases. 

In considering gauge-mediated SUSY breaking (GMSB), we adopt the parameters of 
model I studied for the SUSY/Higgs Run II workshop with tan/? = 2.5, /i > 0, one 
messenger SU(5) generation, and a messenger scale M = 2A. Parameter A is the scale of 
SUSY breaking. We examine six cases in which A varies from 40 to 150 TeV. GMSB does not 
favor associated production at Tevatron energies because it results in a pattern of gaugino 
masses in which M^{M)/M2{M) = as{M)/a2{M), where M3 and M2 are the masses of the 
gluino and weak gaugino, and 0:3 and 0:2 are the SU(3) and SU(2) gauge couplings. Since M 
is a low scale, as is still quite strong. The gluino is generally heavy compared to the other 
gauginos. Selecting A = 40 TeV, Mi, M2, /x = 56.47, 112.8, 241.7 GeV, and we obtain 
neutralino masses m^o ^ = 45, 88, 245, and 281 GeV with m^o < inside a polarization 
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sum. The charErino masses m-± are 82 and 277 GeV and again almost degenerate with 
the masses of X2 4- The gluino mass rrig = 367 GeV, and the squark mass is 471 GeV. The 
spectrum of masses is similar to that of our default SUGRA scenario, and the masses grow 
as we increase A. For comparable gluino and gaugino masses, we find that LO cross sections 
are roughly a factor of 5 (3) smaller at Tevatron (LHC) energies than in the SUGRA model, 
related to the larger squark mass. 



Our anomaly mediated (AMSB) scenario is based on the work in Ref. ||T9|]. It is less well- 
defined in the sense that scalar masses are not understood, and thus the value of fi is not 
determined through radiative electroweak symmetry breaking. However, Mi,M2, and M3 
are well specified. We fix the squark masses to be 350 GeV and choose tan/5 = 2. The 
gaugino masses are controlled by the gravitino mass. The gluino tends to be heavy in this 
scenario, disfavoring associated production. However, the gluino mass has phase vr relative 
to Ml and M2, resulting in constructive interference at LO in the production of Xi; X2) 
and and negative interference in production of x^, in contrast to the SUGRA scenario. 
In AMSB the lightest neutralino is always a W and has a large coupling to (s) quarks, in 
contrast to the B-like lightest neutralino of the SUGRA model. We vary the gravitino mass 
parameter 7713/2 from 30 to 60 TeV. For Mi, M2, and fi = 272, 80, and -300 GeV, we obtain 
neutralino masses 7n^o ^ = 91, 269, 309, and 371 GeV with m^o < inside a polarization 
sum. The chargino masses are m-± = 91 and 318 GeV. The gluino mass rrig = —672 GeV. 
The masses grow as we increase ^3/2- For the gxi channel the LO cross section is roughly a 
factor of 15 larger at Tevatron energies than in the SUGRA model, for comparable masses. 
However, the fact that the combination of the gluino and neutralino masses exceeds 750 
GeV makes this model an unlikely candidate for discovery at the Tevatron. 

Gaugino dominated boundary conditions offer another interesting possibility, exem- 
plified by gaugino-mediated SUSY breaking (gMSB) [l^Jl^ . In this class of models, gauge 
fields propagate freely in the five-dimensional bulk, whereas fermions are confined to one or 
more four-dimensional hyper- surfaces. Supersymmetry is broken at a distant point in the 
extra dimension, giving mass to the gauginos, and the gauge interactions communicate this 
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breaking to the scalar fermions as well. We present results based on the model of Ref. [jlH 



which combines this simple mechanism of SUSY breaking with a model of quark masses and 
mixings. This model has two input parameters, mi/2 and the mass of the down-type Higgs 
at the GUT scale, mn^- Fixing, for example, mn^ = 200 GeV and mi/2 = 150 GeV (which 
determines tan/5 ~ 15 |TB[), we find a spectrum with gluino mass 379 GeV, neutralino 
masses 57, 103, -224, and 249 GeV, chargino masses 101 and 251 GeV, and squark masses of 
about 330 GeV. This spectrum is similar to the default SUGRA scenario, with the principal 
differences that the Higgsino masses are somewhat lighter because the non-universal bound- 
ary conditions at the GUT scale result (through the requirement of radiative EWSB) in a 
non-SUGRA /i-term at the weak scale. For the light gaugino-like states the cross sections 
are virtually unchanged with respect to the SUGRA scenario, whereas the LO cross sections 
for Xg, X4, and xf 5, and 5 times larger than those in the default SUGRA scenario, 

because of the increased phase space for these states. 

An intriguing scenario is that of a light gluino LSP . Gluino masses in the range of 25 
to 35 GeV may still be allowed [E^]. Since the work of Ref. PB[ treats only the strong SUSY 
sector, we make some assumptions about the weak parameters, respecting LEP limits pi 



on neutralino and chargino masses. We choose nig = 30 GeV and rriq = 450 GeV. For the 
weak sector, we adopt masses typical of SUGRA models, discussed above. Since the gluino 
mass is light, there is much more phase space available, and cross sections for associated 
production are substantial at Tevatron energies, reaching ~ Ipb^^ for a wide range of values 
of mi/2. 

Because the GMSB, ^MSB and AMSB cross sections at LO are not too dissimilar from 
those of the SUGRA case at Tevatron energies, we focus our NLO work on the SUGRA and 
light gluino models. 
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III. NEXT-TO-LEADING ORDER CONTRIBUTIONS 



The next-to-leading order contributions to the associated production of gluinos and gaug- 
inos can be separated into virtual corrections that contain internal loops of colored particles, 
and 2-to-3 parton real emission contributions in which a light gluon or quark is emitted. The 
kinematics of the virtual contributions are identical to the Born case described in Sec. |I| 
whereas the presence of an additional out-going particle in the emission contributions re- 
quires integration over a three-body (rather than two-body) phase space. It is useful to 
separate the real emission contributions into parts in which the additional parton's energy 
approaches zero (and thus the three body final state effectively becomes a two-body one) 
and parts in which the additional parton is hard (energetic). We refer to these two parts as 
soft emission and hard emission contributions, respectively. 

All of these NLO contributions contain singularities. The virtual corrections contain 
ultraviolet (UV) singularities that may be absorbed into the definitions of the couplings and 
operators in the usual renormalization procedure. Both virtual and emission contributions 
contain infrared (IR) singularities when the energy of the produced or exchanged particle 
approaches zero. These singularities cancel when the virtual and emission contributions are 
combined. Finally, there are collinear singularities when the produced particle is emitted 
collinearly with another massless colored object. These singularities are absorbed into the 
(universal) NLO definition of the parton distribution functions. 

A. Virtual Corrections 

In this subsection we present the virtual corrections to the associated production of 
gauginos and gluinos in hadron collisions. They arise from the interference of the Born 
matrix elements presented in Sec. |I|with the one-loop amplitudes shown generically in Fig. |^. 
In these diagrams, the crossed regions indicate contributions from self-energy corrections 
(Fig. 1^) and vertex corrections (Fig. ^) that are present one at a time at next-to-leading 
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order. Additional contributions arise from the box diagrams in Fig. ^. We include the full 
supersymmetric spectrum of strongly interacting particles in the virtual loops, i.e. squarks 
and gluinos as well as quarks and gluons. 

Since these virtual loop contributions contain ultraviolet and infrared singularities, we 
regularize the cross sections by computing the phase space and matrix elements in n = 4 — 2e 
dimensions. We then obtain the virtual differential cross section from 



~ t2U2 — 



d% _ nS, 



dt2 du2 Til — e) 



jj? s 



Q{t2 U2 - 1112 s)&{s - {nil + m2) ) 



X(5(s + t + M - - m^) {M^M^* + M^M^*). (21) 

As in the Born case, the matrix elements are summed (averaged) over the colors and spins 
of the outgoing (incoming) particles. 

We calculate the traces of Dirac matrices with the help of the computer algebra program 



FORM [^] using the so-called "naive" 75 scheme. In this scheme, 75 anticommutes with 
all other 7^ matrices, which is justified for anomaly-free one-loop amplitudes |^. The 
75 matrix enters the calculation through both the quark-squark-gluino and quark-squark- 
gaugino Yukawa couplings. The integration over the internal loop momenta is simplified by 
reducing all tensorial integration kernels to expressions that are only scalar functions of the 
loop momentum |^T|. The resulting one-, two-, three-, and some four-point functions were 
computed in the context of other physical processes 0. However, two previously unknown 
divergent four-point functions are computed here for the first time due to the fact that the 
final state particles, i.e. the gluino and the gaugino, have different masses in general. The 
absorptive parts are obtained with Cutkosky cutting rules and the real parts with dispersion 
techniques. The results are collected in Appendix p. 

The virtual one-loop corrections contain ultraviolet divergences that appear as poles in 
1/e in the one- and two-point functions. They are removed by renormalization of the coupling 
constants in the modified-minimal-subtraction scheme (MS) scheme at the renormalization 
scale /i and of the masses of the heavy particles (squarks and gluinos) in the on-shell 



scheme. A difficulty arises from the fact that gluons have n — 2 possible polarizations, 
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whereas gluinos have 2, leading to violation of supersymmetry in the MS scheme. The 
simplest procedure to restore supersymmetry, which we adopt here, is through a finite shift 
in the quark-squark-gluino Yukawa coupling: 



9s = 9s 



as 



1 + — I -A^ 
47r V 3 



1 



9s 



1 + 



as 



SvrJ 



(22) 



This shift was discussed first in Ref. ||33 |. 

The virtual corrections can be classified into aCp and a Nc color class depending on the 
color fiow or the Abelian or non-Abelian nature of the correction vertices. In addition to UV 
singularities they have collinear and infrared singularities that appear as 1/e or poles in 
the derivatives of the two-point- and in the three- and four-point functions. The generally 
divergent scalar integrals are always multiplied by finite coefficient functions proportional 
to parts of the Born matrix elements. The full result is given in Appendix 



B. Real Emission Contributions 

At NLO, the production of gluinos and gauginos receives contributions from real emission 
of gluons or massless quarks and anti-quarks. In the following sub-sections both of these 
types of two-to-three partonic contributions are dealt with separately. 



Following the notation developed in Ref. |24[, we express our results in terms of the 
following sets of invariants. 



S = {Pa+Pbf 

S3 = {P3+P2? - ml 

t= {Pb- P2f 
U = {Pa -P2Y 

ue = {Pb - Pif - ml 



Sb = (Pl +P2f 

Si = {Ps+PiY - ml 
t' = {Pb - PsY 

U' = {Pa-p-if 

uj = {Pa-Pif - ml, 



(23) 



where pa and ph are the four-momenta of the incoming (massless) partons, pi and p2 are the 
g and x momenta, and p^ is the momentum of the additional massless parton. We also find 
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it useful to define tlie following derived quantities: 



tl 


= t — ml 


t2 


= t — m\ 


Ml 


= u — m\ 


U2 


= u — ml 


A„ 


2 2 

= m, — m= 


At 


= m\- m? 




= Uq + Au 




= M7 + At 


S4A 


= S4 + Au 


S3A 


= S3 + At . 



(24) 



Energy and momentum conservation provide relations among these quantities: 

54 = S + ^2 + ^il S3 = S + U6 + ^7 (25) 

55 = s + t' + Me = — s — ^2 — ^' 
= —s — U2 — u . 

The n-dimensional three-body phase space may be derived conveniently if we integrate 
the general fully differential cross section in the 1-3 rest frame |3^. In this frame the 



4-dimensional components of the n-dimensional momenta are expressed as: 

Pa = (Wa, 0, LOa siu U)^ COS -0) (26) 

P6 = (wft, ,0,0,a;f,) 

Pi = (£^1, — ciJs sin dx sin Q^-, —'^z sin cos 6^2, — a'3 cos ^1) 

P2 = (-E2,a;aSin?/',u;aCos?/' + uj^ 

P3 = (a;3 , uj-i sin 6^1 sin 6*2 , ^73 sin 6^1 cos 6^2 , 0^3 cos 6'i ) , 



with 



S + U2 s + t2 

Ua = , ^ UJb = — , ^ (27) 

2ys4 + mf 2A/s4 + mf 

S4 S4 + 2 

1^3 = — / ^1 



2a/s4 + 2a/s4 + 

^ t2 + M2 + 2m? -s(s4 + m? + mg) + t2^i2 

i^2 = , COS-W = -. TT : 

2^54 + ^2 (s + t2)(s + U2) 
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Using this parameterization, we may express the invariants defined in Eq. (^) in terms of 6i, 
02, and the 6'-independent variables uj(^a,b,3), -^(1,2)) i^- the real emission contributions, 
it is sometimes convenient to parameterize these momenta with the z axis aligned along pa 
or P2. As these alternate frames are related by a simple spatial rotation, the expressions for 
£'(12), and uj(a,b,3} remain unchanged. The general three-body cross section may be expressed 
in this frame as 



t2U2 — S 1712 



S fl"^ 



e{t2U2- smi,)&{s4) (28) 



ds4^du2dt2 2s'^T{l — 2e) 

r 2 

X e(s - (mi + 1712^) 7 ^^r- 5{S + t2 + Ui- S4) / dVlnlM^] , 

(s4 + mf)^-^ J 



in which \A4^\ is the real emission matrix element squared, summed over final spins and 
colors and averaged over initial spins and colors, and the ra-dimensional angular integration 
is c/fi„ = sin^-2e(^^^^) sin-2^(^2) ^^^2- 

In evaluating the integration over the angular variables in Eq. (pH]), we follow the proce- 
dure outlined in Ref. |3^, in which we use the relations among the invariants, Eq. (pSf), to 



reduce all of the angular integrals to the form. 



r sm''''i9^)d9, r sm-^'i92)d92 (29) 
Jo Jo 

x{a + bcos9iy'{A + Bcos9i + Csm9i cos6'2)"'. 



the necessary expressions for which may be found in Ref. [^. The angular integrations 
involving negative powers of t' and u' produce poles in e which correspond to collinear 
singularities in which particle 3 is collinear with particle a or 6 (c.f. Figs. (^ and (|^). 
Because these singularities follow a universal structure, they may be removed from the cross 
section and absorbed into the parton distribution functions according to the usual mass 
factorization procedure ||3^ . The non-zero mass of the gluino kinematically forbids collinear 
emission, and thus the gluino has no associated collinear singularities. 
The collinear singular pieces have the factorized form 
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dt2 dU2 



I dxiXi I dx2X2 rfcj(xi, 
-^0 "^0 k,l 



(30) 



(it2 du2 

in which s = Xi X2 s, U2 = xi U2, ^2 = X2 ^2- The universal sphtting functions Tij contain the 
colhnear divergences associated with incoming parton j sphtting into parton i, and the hard 
scattering cross section, d'^d^i/di2 du2, is free from singularities. The splitting functions may 
be redefined by an arbitrary finite term, and thus one must choose a factorization scheme. In 
order to use recent sets of parton distributions extracted from data we adopt the MS scheme, 
in which the splitting functions at 0{as) are 



rij(x,/i|,/i^e) = 6ij6{x- 1) + ^ 



-- + 7E-log(47r)+logf^ 



P,,{x) 



(31) 



In the above expression, the Pij{x) are the Altarelli-Parisi evolution kernels [p6 

Pqq(yXi) = Pqq{Xi) = Cp 



l±^e(l - 5, - Xi) + (2 log 5, + \]5{l- Xi) 



(32) 



Pgqi^Xi) 



Pgqi^i) 



l + il-Xi 



Xi 



Pqg{Xi) = Pqg{Xi) = Tp [x^ + (1 - Xif] 



Pgg{Xi] 



2Nc 
+ 



Xi { 1 Xi , 



+ Xiil - Xi) -2 



Xi - 6i 



2Nclog6, + =;P^ 



1 

2' 



6{1 - Xi 



with Cf = 4/3 and Tp = 1/2 for Nc = 3 colors of quarks; (3^ = llNc/Q - 2nFTF/S. 
The quantities 6i express the slicing of S4 into hard and soft regimes in terms of the Xa 
and Xb variables. They can be related to the A of Section |IIIB 1| hj 6a = A/ {s + U2) and 
6b = A/(s + t2). 

We set the renormalization and factorization scales equal to each other, fip = IJ', and 
expand Eq. ([30| ) to 0{as) to derive the expression for the reduced cross section. 
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dt2 du2 dt2 dU2 

d'^a'^j{xis,t2, X1U2, n"^) 



(33) 



+ / dxiXiPii{xi) 

27r e Jq 



dt2 du2 



+ t dX2X2P,,{x2) d'-U^-'^^^t2,U2,f.^) 



2vr e Jq dt2du2 
We employ the compact notation = — 7_e + log(47r); d^a^^./ dt2du2 is the leading 
order cross section for ik gx- The resulting hard scattering cross section is free from 
collinear singularities, as the implicit e-dependence of d^aij/dt2 du2 cancels with the explicit 
e-dependence of the second and third terms. 



1. Gluon Emission 

The NLO real contributions with an additional gluon in the final state, 

^ 9 9X, (34) 

proceed from the Feynman diagrams shown in Fig. ^. As was the case for the leading 
order cross section for production of neutralinos with gluinos, the set of graphs in which a 
right-handed squark is exchanged cannot interfere with the graphs in which a left-handed 
squark is exchanged because the incoming quark and anti-quark must have definite helicity. 
Production of charginos in association with gluinos involves only left-handed squarks. 

In addition to the collinear singularities described above, this set of corrections also has 
infrared singularities that arise when the energy of the emitted gluon approaches zero. These 
singularities appear as poles in S4 in the reduced cross section, and must also be extracted 
so that they can be combined with corresponding terms in the virtual corrections and shown 
to cancel. 

To make this cancellation conveniently, we slice the gluon emission corrections into hard 
and soft pieces, 

~ c^g4 , , , + / , , , > (35) 



dt2 dU2 Jq dt2 dU2 ds4 J/^ dt2 dU2 dS4, 
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where A is an arbitrary cut-off between what we call soft gluon radiation and hard gluon 
radiation. When the cut-off is much smaller than the other invariants, A ^ s, t, u, m|, the 
S4 integration for the soft term becomes very simple and can be evaluated analytically. This 
operation results in singular terms 



2tt J \ \ smf J J e 

In Eq. (|36|) we distinguish the contributions from the Nc and Cp color classes. 

The singular expression may then be combined with the virtual corrections discussed in 
Sec. [Ill A| to yield the combined "soft and virtual" contribution free from infrared singular- 



ities. The residual finite soft contributions are presented in Appendix 0. 

In the hard gluon regime, there are collinear singularities, but no IR singularities, and 
thus the most singular terms are proportional to e^^. After the mass subtraction described 
above is performed, the results are singularity- free, and they can be presented as minimally 
subtracted, singularity-free integrals, 

Hf (01,02)) = r de, rrf^2sini-2^(^i)sin-2^(^2)/(^i,^2) (37) 
Jo Jo 

-i Jim (y^ dO^ d92sm^-''{e,) sin~2^(^2) /(^i, ^2)) , 

which contain only the finite terms with the poles subtracted. The resulting expression 
consists of a simple power series in e, which may then be evaluated in 4 dimensions by setting 
e 0. Note that the function f{6i,62) can involve coefficients for angular expressions that 
have mass dimension, and thus the mass dimension of I{f{6i,62)) will depend on f{6i,62)- 
The gluon emission matrix elements are presented in Appendix 0. 

The cutoff on the S4 integration introduces an implicit logarithmic dependence on A that 
is matched by the explicit logarithms of A which appear in the combined soft and virtual 
term. The total correction is independent of the value of A. Choosing for illustration 
^1/2 = 400 GeV, we display in Fig. ^ the dependence of various contributions on our 
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cutoff S. The Bom contribution is obviously independent of S, but its contribution helps to 
show the relative magnitude of different terms. The combined soft and virtual contribution 
is positive but falls as an explicit analytic function of log5. The hard part of the gluon 
emission contribution is negative, but its numerical value grows more positive as a implicit 
function of \og6. The figure shows that two contributions balance each other well, and the 
combined soft/virtual plus hard contribution is independent of the cutoff for 6 < 2 10^^. 
The figure also shows that the net small next-to-leading order contribution is obtained after 
large cancellations take place. The case chosen for display in Fig. ^ is a worst case. With 
= 400 GeV, m^o = 679 GeV, and rrig = 1012 GeV. The energy chosen is that of the 
Tevatron, = 2 TeV, so phase space limitations are relatively severe for this set of masses. 
For all other cases, the cross sections are independent of 6 for a larger range of 6. 

2. Light Quark (Anti- Quark) Emission 

A second set of real emission corrections involves an additional light quark (or anti-quark) 
in the final state, through the partonic reactions, 

(19 ^ <19X,Q9 ^ Q9X- (38) 

The set of Feynman diagrams contributing to emission of a quark is shown in Fig. 0. They 
include diagrams in which an incoming gluon splits into a qq pair as well as diagrams in 
which an intermediate squark splits into a quark and either a gluino or gaugino. This 
set of corrections does not have an IR divergence, and thus it is not necessary to slice it 
into hard and soft regimes. However, after all of the initial state coUinear singularities are 
removed by the mass factorization procedure described above, the matrix elements may 
still contain integrable singularities if the mass of the squark is larger than the mass of the 
gluino or gaugino. In these cases, the intermediate squark state can be on its mass-shell, 
and the variables s^a and s^a go to zero inside the region of integration. This problem was 
encountered previously [0, and we follow the same procedure. These singularities represent 
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the LO production of a squark and a gluino or gaugino, followed by the LO decay of the 
squark. They may be removed if one includes the full Breit-Wigner form for the squark 
propagator, which regulates the squark resonance by the squark width. This procedure 
amounts to the replacements, 

— -^—J^ ^ ^ (—) - i ^ (39) 



where V indicates the principal value function, and the final distribution identity holds in 
the limit of small squark widths, Vq <^ rUq. The replacement removes the singularities, and 
when both S4a and S3A are zero it generates an additional real term from the product of 5 
functions. 

There is a further subtlety associated with the requirement that we not double-count the 
region of phase space in which the squark is on-shell. Properly, the kinematic configuration 
with an on-shell squark is included in the LO production of a squark and a gluino or a squark 
and a gaugino, and thus should not be considered as a genuine higher order correction to 
the production of gluinos with gauginos. To avoid double counting, we thus subtract the 
on-shell squark contribution by defining the total cross section (for illustration, we deal with 
the S4A singular case). 



dSi I dt2 - — ; — = / ds4 — 5 — 5-. (40) 

Jt-{si) dtidsA Jo s|A + "ig~r-~ 

The on-shell contribution then corresponds to /(O) / (s|a + ^|r|), with 



/(O) = 4, !Vi £lzlM ^ BR{q q~g) {sl^ + m? T?) 5(.4a). (41) 

TT in 



9 

It can be subtracted leaving a genuine NLO contribution 

.NLO_ r^'^\. /Csia)-./(0) 



a 



r -/g) (42) 



which may once again be expressed as a principal value function, since S4A / (s4a + ^ti? F?) — >■ 
'P(1/'S4a) in the limit of small squark width. The ssa singular terms may be treated in a 
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similar way, with the added comphcation that the integration over s^a is hidden in the 
angular integrations [j^. 

The quark emission matrix elements are presented in Appendix 0. 

For the parameters of the SUGRA model that we adopt, the gluino mass remains greater 
than the squark mass for all values of mi/2, and there is never an intermediate squark-to- 
gluino-plus-quark final state singularity. However, the two chargino masses and all four 
neutralino masses are always less than the squark mass, and the final-state on-shell squark- 
to-gaugino-plus-quark singularity comes into play in all cases. In the light gluino model, 
with rriq = 450 GeV, the gluino mass of 30 GeV is light enough that on-shell intermediate 
squark decay into a gluino is always active. In this light gluino model, as mi/2 is varied from 
100 to 400 GeV, the masses of the lighter gauginos (x?, xt^ and x^) remain less than the 
squark mass so that the on-shell intermediate squark decay into a gaugino is active over the 
whole range of mi/2 for these light channels. The situation changes for the heavier gauginos 
(Xg, x1, and xf)- For small mi/2 their masses are below the squark mass, and the on-shell 
decay is active. However, above roughly mi/2 = 250 GeV, the masses of the heavier gauginos 
exceed the squark mass and the on-shell possibility closes. 



IV. QUANTITATIVE RESULTS 

In this section we collect our main results on total and differential cross sections for the 
associated production of gauginos and gluinos at Tevatron and LHC energies. 



A. Scaling Functions 

We begin with the cross section at the parton level expressed as 

= + 4vr«5(/i) Ur^-^'iv, /i) + fSiv, /i)] } • (43) 

It has been integrated over the Mandelstam invariants t and S4 and depends on the partonic 
center-of-mass energy s through the scaling variable 
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4m^ 

where m is the average mass of the produced sparticles 

mi + 7712 

m = . 



' -1, (44) 



(45) 



It also depends on the produced masses mi and m2 and on the squark mass rriq (through 
the internal squark propagator). The common renormalization and factorization scale is 
denoted by fi. The partonic initial state is labeled i,j = g, q, q. 

Equation ( P3| ) defines the dimensionless scaling functions fij, studied in Ref. ^7j. These 



functions are independent of the coupling constants a and as, of parton densities, and of 
the collider type and energy. They permit precise checks of individual contributions and of 
the threshold, resonance, and high energy behaviors of the production process. 

The Born and the summed virtual and soft scaling functions fij~^^ receive contribu- 
tions only from qq initial states, where q = u or d, with the possible emission of a soft and/or 
collinear gluon. The hard scaling function has contributions from qg initial states when 
an additional quark or antiquark is emitted together with the gluino and the gaugino. We 
eliminate the explicit dependence of the soft contributions on the technical cut-off 6 = A/m^ 
by subtracting the log^ ' ' 6 terms. These terms are then added to the hard contribution 
such that this contribution is also independent of 6. In Sec. |T| we show that our results are 
independent of 6 at least in the range 6 G [10^^; 10"'^], and we use the value 6 = 10~^ in the 
following. 

The scaling functions for the production of a ^ and a X2 presented in Fig. ^ and those 
for a g and a xf in Fig- HO- Here we set the scale /i equal to the average particle mass m. The 



masses are those of our default SUGRA scenario. As discussed in Sec. [Ill B 2[ the emission 
of an additional quark or antiquark can lead to intermediate on-shell squarks and therefore 
to a singular squark propagator in Feynman diagrams. After the LO two-body q + g g + q 
contribution is removed, the remaining integrable singularities can be identified as spikes 
in the gu and gd scaling functions in the two figures. They yield finite contributions after 
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integration over the momentum fractions Xa and Xb of the incoming partons or, equivalently, 
over the partonic center-of-mass energy s = XaX^S. 



Evident from Figs. |^ and [1^ is that next-to-leading order contributions do not aher ei- 
ther the threshold or high energy asymptotic behaviors in rj, unlike, e.g., the situation for 
pair production of heavy quarks [^,^. The combined virtual and soft scaling functions 
contribute negatively but are small in magnitude when compared with the hard scal- 
ing functions f^. The figures show that one should expect only modest enhancements in 
predicted rates when the uu, dd, and ud production channels are dominant, as is true at 
Tevatron energies, for which the quark/ ant iquark parton luminosity is large and the range 
in 7] is limited ( 77 < 15 for these two channels in our default scenario). The contribution 
of the qg channel can become important if phase space is open to large values of rj. At 
LHC energies, 77 extends to nearly 800. This range in rj, along with the large qg luminosity, 
suggests that the qg channel will supply significant enhancements in the predicted rates at 
LHC energies, as is demonstrated below. 

Scaling functions for the gxl, gx^, and gxf channels show behavior similar to that seen 
in Figs. 1^ and with the notable exception that the positive excursion at large 77 in the 
qg channels is relatively more prominent for gxl and gxf than for gx2 and gxf- The 
channel is distinguished primarily by the fact that the peak in the rj distribution at the Born 
level occurs near rj = 0.4 whereas the peaks occur at larger rj, in the range of 77 = 1 to 2, in 
all other cases. There is also a noticeable difference in threshold behaviors of the Born and 
NLO hard-gluon emission contributions for this channel. We relate these differences to the 
fact that only the gxs channel exhibits positive interference at the Born level between the 
t- and u-channel contributions, c.f., Eq. (^). Because tti^o is negative in the SUGRA model, 
the term proportional to Xf^ in Eq. is positive for gx^ production but negative in all 
other cases. 
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B. Hadronic Total Cross Sections 



The hadronic total cross section is obtained from the partonic cross section through 

a^'^^{S,^)= f dXa f dXbf!''{Xa,fi)fl^^{Xb,fi)aij{XaXbS,fi), 



where 



(46) 



r='4, (47) 



and a/S* is the hadronic center-of-mass energy (2 TeV for Run II at the Fermilab pp collider 
Tevatron and 14 TeV at the CERN pp collider LHC). Our NLO predictions are calculated in 
the MS scheme with the CTEQ5M parametrization for the parton densities f{x,fi) in 



the proton and antiproton and a two-loop approximation for the strong coupling constant 0:5 
with A*^^^ = 226 MeV. This value of A^^-* is used also in the renormalization group evolution 
equations for our SUSY scenarios. In obtaining LO cross sections, we use the CTEQ5L LO 
parton densities and the one-loop approximation for as, with A^^^ = 146 MeV. 



1. SUGRA model results 



For the SUGRA scenario, we present the total hadronic cross sections for the associated 



production of gluinos and gauginos at Run II of the Tevatron in Fig. |Tl| and for the LHC 
in Fig. |T2|. We vary the SUGRA parameter mi/2 from 100 to 400 GeV and keep the other 
SUGRA parameters fixed at the values described in Sec. P A| . The squark mass runs from 



250 GeV to 890 GeV in this region. The cross sections are presented as a function of the 
physical gluino mass rrig. The corresponding gaugino mass ranges are 31 to 163 GeV for xl 
62 to 317 GeV for and xf, 211 to 666 GeV for X3, and 240 to 679 GeV for xl and xf. The 
chargino cross sections are summed over positive and negative charges. The renormalization 
and factorization scale fi is set equal to the average particle mass m. We truncate Fig. ^ 
at a cross section of 10~^pb since the anticipated integrated luminosity at Run II is at most 
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30fb"\ For the convenience of the reader, we provide numerical values of the cross sections 
in Table I for a few selected points. 

For small rrig one might expect the largest cross section for the lightest gaugino, Xi- 
However, its coupling is dominantly of type B and therefore smaller than the W^-tjpe cou- 
pling of X2 which, in turn, has a larger cross section despite its larger mass. The heavier 
gauginos X34 and xf are dominantly Higgsino and are therefore suppressed by several or- 
ders of magnitude with respect to the lighter gauginos because of the light quark Yukawa 
couplings. At the LHC, the xf cross section is dominant. At small values of rrig, the LHC 
cross sections are a factor of about 30 greater than at the Tevatron, and at large rrig, the 
factor is about 10^. 

Comparing the NLO predictions in Figs. |ll|and ^ (solid curves) with the LO predictions 
(dashed curves) we observe that the NLO corrections are all positive and substantially larger 
at the LHC than at the Tevatron. At the Tevatron, some of the NLO predictions fall below 
the LO predictions at large mass, a point to which we return below. The NLO enhancements 
are more evident in the ratio of the NLO over the LO cross section {K factors) shown in 



Figs. |I3| and |T4[ The K factors are computed at the scale fx = m. At the Tevatron, the NLO 
corrections amount to at most a 10% increase in cross section. At the LHC, they appear 
generally in the range of 20 to 40% but can amount to a factor of 2 for Xi and xf. 

The very modest size of the NLO enhancement at the Tevatron is somewhat expected 
from the behavior of the scaling functions, but it is also attributable partly to differ- 
ences in the NLO and LO parton densities. Recalculating the K factors with the CTEQ4 
parametrization |^0|, we find increases in K by as much as 0.1 at the Tevatron energy. The 
change from CTEQ4 to CTEQ5 is interesting. The u quark density at NLO decreased by 
1 to 5% and the d quark density at NLO increased by up to 10% over the range 0.1 < x < 
0.6. On the other hand, the LO u quark density increased by about 1% and the LO d quark 
density by up to 20%. These purely parton density effects result in a net increase in the LO 
cross sections and decrease in the NLO cross sections, a drop in the calculated K factors 
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from CTEQ4 to CTEQ5. 

The contribution from the qg initial state at the energy of the Tevatron is insignificant 
(less than 10^'^ of the total) for all six gaugino channels, but it is considerable at the energy 
of the LHC In Fig. ^ we show the fraction of the NLO cross section at the LHC attributed 
to the qg initial state. At the energy of the LHC, the contribution from the qq channel is less 
dominant, and the qg contribution becomes significant owing to the large gluon density. A 
similar effect is seen in single top quark production for comparable values of the average 
produced mass m. 

The large K factor for x\ ^"^^ X2 production at the LHC is associated with a large 
contribution from the gq channel at this energy. For these two channels, there is strong 
interference at LO between the t- and u-channel exchange diagrams that does not occur in 
the NLO quark emission graphs. 

The dependence of the predicted cross sections on the renormalization and factorization 



scale is reduced considerably at next-to-leading order. As an example, in Fig. |T6| we show the 
scale dependences for the X2 channel at the Tevatron. All six channels show similar behavior 
at the Tevatron. Cross sections vary by ±23% at LO but only by ±8% at NLO as the scale 
ratio /i/m is varied over the range 0.5 to 2.0, a substantial improvement in reliability. Here 
m is the average mass for the default scenario. At the Tevatron, the NLO and LO cross 
sections intersect at scale ratio near unity. In Fig. |17| we show the /i dependences for the X2 
channel at the LHC. Again, all six channels show similar behavior. Cross sections vary by 
±12% at LO but only by ±4.5% at NLO in the region 0.5 < ^/m < 2.0. However, unlike 
the situation at the Tevatron, owing to the important contribution of the qg channel at the 
energy of the LHC, the NLO and LO cross sections do not intersect at scale ratio near unity. 
Instead, if they intersect at all, the crossing point is at a very low scale. 

Uncertainties in the cross section from parton density variation may be estimated roughly 
if we compare NLO results obtained with CTEQ4M and CTEQ5M. For X2 production at the 
Tevatron in our default SUGRA scenario, we compute cross sections of 0.0244 pb (CTEQ4M) 
and 0.0219 pb (CTEQ5M), a 12% difference. At the LHC, the cross sections are 1.138 pb 
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(CTEQ4M), and 1.095 pb (CTEQ5M), a 4% difference. 



2. Light gluino model 

For the light gluino model, we present the total hadronic cross sections for the associated 



production of gluinos and gauginos at Run II of the Tevatron in Fig. |18| and for the LHC in 
Fig. As mentioned in Sec. II A, we fix rrig = 30 GeV and rriq = 450 GeV. We display the 
cross sections as functions of the mass of the common GUT-scale fermion mass mi/2- As we 
vary mi/2 from 100 to 400 GeV, the gaugino masses range from 31 to 163 GeV for Xi^ 62 
to 317 GeV for and xf, 211 to 666 GeV for xl and 240 to 679 GeV for xl and xf. It is 
worth noting that the coupling strengths also vary with mi/2- The chargino cross sections 
are summed over positive and negative charges. The renormalization and factorization scale 
H is set equal to the average of the masses of the gaugino and gluino in each of the channels. 



Evident in Figs. |T8| and |T9| is that the cross sections do not depend strongly on the 
gaugino masses. The values of mi/2 in these figures extend below the value ^150 GeV 
believed excluded since LEP data |^T| set a lower bound on the mass of x^ of about 100 
GeV. Nevertheless, even above mi/2 = 150 GeV, the Tevatron cross sections for the three 
lighter gluino channels are predicted to be in the range of 0.1 to 0.5 pb. The cross sections 
would be increased if rriq were reduced from the value 450 GeV that we use. 

The relatively large cross sections suggest that associated production is a good channel 
for discovery of a light gluino at the Tevatron, for closing the window on this possibility, 
and/or for setting limits on light gaugino masses. We remark in this connection that the 
usual searches for a light gluino LSP begin with the assumption of pair production of gluinos. 
In this situation, the dominant background is QCD pair production of hadronic jets. Hard 
cuts on transverse momentum must be made to reduce this background to tolerable levels. 
The cuts, in turn, mitigate against gluinos of modest mass. By contrast, if light gluinos are 
produced in association with gauginos, one can search for light gluino monojets accompanied 
by leptons and/or missing transverse energy from gaugino decays. 



29 



The K factors for the hght gluino case are shown in Figs. ^ and ^ At the Tevatron, 
the NLO corrections amount to a 30% to 40% increase in cross section. At the LHC, they 
are generally in the range of a factor of 2 to 3.5. The large K factors owe their origins to 
the important role of the gq channel. The gluon parton density is very large at small values 
of X. Contributions from the gq production channel are more intense in the light gluino case 
than in the SUGRA case where average produced masses and, thus, typical values of x are 
larger. At the energy of the Tevatron, the gq channel accounts for more than 20% of the 
NLO cross section in the Xi and Xi channels at small values of mi/2 and more than 10% 
at large mi/2- At the energy of the LHC, the fraction exceeds 60% in the Xi-, X21 and xf 
channels for the entire range of mi/2- It hovers above 50% in the Xa; x\i and channels 
at small values of mi/2 and near 20% at large mi/2- 

Motivated by the curious behavior of the K factors for the three heavier gaugino channels 
at the LHC energy, we examined the change in the scaling functions for the x\ case as mi/2 is 
varied. In the qq channel, the net (virtual plus soft plus hard) NLO contribution is positive 
for all ?7, and, relative to the Born contribution, its magnitude grows gradually with mi/2- 
The qq channel accounts for a slightly increasing component of K, hovering about 1.5. On 
the other hand, the contribution from the gq channel changes markedly as mi/2 increases. 
Below about mi/2 = 200 GeV, its scaling function is large, with significant support below 
7] = 1, where the gluon parton density is large. As mi/2 increases above 200 GeV, the gq 
scaling function decreases in magnitude. The gq channel supplies a component of K that 
increases slightly from about 2 to 2.5 as mi/2 increases to 200 GeV and then falls gradually 
to below 0.5 at mi/2 = 400 GeV. 

We attribute the sharp decrease of the K factor for the three heavier gaugino channnels 
to the role of on-shell intermediate squark decay into a gaugino, discussed in Sec. [IIIB 2 



In the light gluino model, as mi/2 is varied from 100 to 400 GeV, the masses of the lighter 
gauginos (xi, X2) and xf ) remain less than the squark mass so that the on-shell intermediate 
squark decay into a gaugino is active over the whole range of mi/2 for these light channels. 
The situation changes for the heavier gauginos (xs, Xai and X2')- small mi/2 their 
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masses are below the squark mass, and the on-shell decay is active. However, above roughly 
mi/2 = 250 GeV, the masses of the heavier gauginos exceed the squark mass and the on-shell 
possibility closes. 



C. Differential Cross Sections 



In Figs. ^ and ^ we display the differential cross sections in the rapidity y and in the 
transverse momentum pt of Xi at the energy of the Tevatron collider. Here the xf and g 
masses are set at their default values in the SUGRA scenario, 101 and 410 GeV, respectively. 
The rapidity distribution is obtained after integration over all pt-, and the pt distribution 
after integration over all y. More restrictive selections could be made, but until experimental 
conditions are better known, any such restrictions (cuts) would be unmotivated. The NLO 
(solid) curve in Fig. ^ shows a modest enhancement in the rapidity distribution in the 
central region with respect to the LO (dashed) curve, but the shape of the distribution is 



unchanged qualitatively. The pr distributions in Fig. |23| show that the NLO contribution 
tends to shift the distribution to somewhat smaller values of pt- Since the contribution of 
the gq initial state is very small at the energy of the Tevatron, the shift in the px distribution 
is associated with next-to-leading order corrections in the dominant qq initial state. 

The features of the y distributions are qualitatively similar for all gauginos except for 
the expected and systematic narrowing of the y distribution with increasing gaugino mass. 
We show one example representative of the full set. The pt distributions for the different 
gauginos are also qualitatively similar except that the maximum in the distribution moves 
to larger pt as the gaugino mass is increased. The location of the peak is specified roughly 
by m/2, where m is the average mass of the produced sparticles. The one exception is Xs 
production. In this case, the peak occurs at a smaller value (about 100 GeV for the default 
masses), an effect correlated with the fact that the location of the peak in the scaling function 
occurs at a smaller value of r]. Interference effects enhance the cross section at small pt for 
x\ production. 
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We show differential cross sections in y and pt for Xi tlie energy of the LHC coUider 
in Figs. ^ and ^ and in Fig. ^ we present the pr distribution for x^. The rapidity 
and transverse momentum distributions are much broader than at the Tevatron. As at 
the energy of the Tevatron, the features of the y distributions are quahtatively similar for 
all gauginos except for the expected and systematic narrowing of the y distribution with 
increasing gaugino mass. The NLO contributions enhance the y distributions at all y. 

At the energy of the LHC, the px spectra are qualitatively similar for the relatively light 
Xi, X25 ^iid Xi states, illustrated by the Xi case in Fig. The enhancement factor K is 
near unity at small pt but becomes sizeable at larger pt- For x\ ^-nd X2 , the px spectra are 



altered significantly by NLO contributions. As illustrated in Fig. the NLO contribution 
associated with the important qg channel fills in the distribution at small px and softens 
the overall px distribution in these two cases. At LHC energies, it is therefore not a good 
approximation to assume that the enhancement factor K is roughly independent of px- 



V. CONCLUSIONS 

In this paper we report a complete next-to-leading order analysis of the associated pro- 
duction of gauginos and gluinos at hadron colliders. If supersymmetry exists at the elec- 
troweak scale, the cross section for this process is expected to be observable at the Fermilab 
Tevatron and/or the CERN LHC. It is enhanced by the large color charge of the gluino 
and the relatively small mass of the light gauginos in many SUSY models. Associated pro- 
duction represents a chance to study in detail the parameters of the soft SUSY-breaking 
Lagrangian. The rates are proportional to the phases of the gaugino and gluino masses, 
and to the mixings in the squark and chargino/neutralino sectors. Thus, in combination 
with other channels, associated production could allow one to measure some or all of these 
quantities. 

The physical gluino and gaugino masses that we use, as well as the gaugino mixing 
matrices, are based on four popular SUSY breaking models plus a fifth scenario in which the 
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gluino mass is relatively light. Because the LO cross sections in gauge-mediated, gaugino- 
mediated, and anomaly-mediated supersymmetry breaking models are not too dissimilar 
from those of the SUGRA case at Tevatron energies, we focus our NLO work on the SUGRA 
model and on a model with a hght gluino LSP, with rrig — 30 GeV. 

In the SUGRA model, the largest cross sections at the Fermilab Tevatron energy are 
those for X2, enhanced by its H^a-like coupling with respect to the 5- like Xi^ the xf, 
which is about equal in mass with the xt- The NLO corrections to associated production 
are generally positive, but they can be modest in size, ranging in the SUGRA model from a 
few percent at the energy of the Tevatron to 100% at the energy of the LHC, depending on 
the sparticle masses. In the hght-gluino case, NLO contributions increase the cross section 
by factors of 1.3 to 1.4 at the energy of the Tevatron and by factors of 2 to 3.5 at the energy 
of the LHC. The large K factors owe their origins to the important role of the gq channel 
that enters first at NLO. 

Owing to the NLO enhancements, coUider searches for signatures of associated produc- 
tion will generally discover or exclude sparticles with masses larger than one would estimate 
based on LO production rates alone. More significant from the viewpoint of reliability, the 
renormalization and factorization scale dependence of the cross sections is reduced by a 
factor of more than two when NLO contributions are included. 

At Run II of the Fermilab Tevatron, for an integrated luminosity of 2 fb~^, we expect 
that 10 or more events could be produced in each of the hghter gaugino channels of the 
SUGRA model, ^Xi? 9X2^ ^^(^ dxf^ provided that the gluino mass rrig is less than 450 GeV. 
The cross sections for the three heavier gaugino channels, gx^, gx% and gxf, are smaller 
by an order of magnitude or more than those of the lighter gaugino channels. In the light 
gluino model, more than 100 events could be produced in the three lighter gaugino channels 
provided that the common GUT-scale fermion mass mi/2 is less than 400 GeV, and as many 
as 10 events in the three heavier gaugino channels as long as mi/2 is less than 200 GeV. At 
the higher energy and luminosity of the LHC, at least a few events should be produced in 
every channel in the SUGRA model and many more in the light gluino model. 
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The shapes of the rapidity distributions of the gauginos are not altered appreciably by 
NLO contributions, but the locations of the maximum cross section in transverse momen- 
tum (pt) are shifted to smaller values by NLO contributions. At LHC energies where the 
contribution of the qg initial state is important, modifications of the pr spectra can be 
pronounced. 

The relatively large cross sections suggest that associated production is a good channel 
for discovery of a light gluino at the Tevatron, for closing the window on this possibility, 
and/or for setting hmits on hght gaugino masses. The usual searches for a light gluino 
LSP are based on the assumption that gluinos are produced in pairs. In this situation, 
the dominant background is QCD production of hadronic jets. Hard cuts on transverse 
momentum must be made to reduce this background to tolerable levels. The cuts, in turn, 
mitigate against gluinos of modest mass. By contrast, if light gluinos are produced in 
association with gauginos, one can search for light gluino monojets accompanied by leptons 
and/or missing transverse energy from gaugino decays. 
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APPENDIX A: VIRTUAL LOOP CONTRIBUTIONS 



At the one-loop level in SUSY-QCD, virtual corrections contribute to the hadroproduc- 
tion of supersymmetric particles through the interference of self-energy corrections, vertex 
corrections, and box diagrams with the tree-level diagrams. For the associated production 
of gluinos and gauginos one has to calculate the self-energy corrections in Fig. ^, the vertex 
corrections in Fig. ^, and the box diagrams in Fig. ^. 

The self-energy corrections for the external quark and antiquark factorize the complete 
Born matrix element and are independent of the underlying scattering process. Since the 
quark and antiquark are treated as massless, only the squark-gluino loop contributes to 



where the factor of two accounts for the sum of the quark and antiquark contributions. 
Functions Bq{p'^; 1711,1712) and i?Q(p^; mi, 7712) stand for the scalar two-point integral and its 
derivative. They are defined in Appendix 0. The quark-gluon loop contains an infrared 
singularity 



where = (47r)'^/(167r^)e~'^'^^. This infrared singularity is not shown in the result above 
since it is canceled by an ultraviolet singularity when infrared and ultraviolet singularities 
are not distinguished in dimensional regularization. 

The external gluino self-energy also factorizes the complete Born matrix element and is 
independent of the underlying scattering process: 



\M'i\ = i[Bo{0- nig, nig) + S^,(0; nig, mg){nil - nil)]2CFgl\M 



(Al) 



(A2) 




= z 



+ i 



+ i 



'"^t) ^ + Bf){m-;nig,nit) 



+ B'oiml; nig, 0)(m| - m|) njih^.lM 



(A3) 
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The scalar one-point integral Ao(m) is defined in Appendix ^ Except for the gluino mass 
and the A'",^ color factor, the gluon-gluino loop is identical to the heavy-quark self-energy. It 
contains an infrared singularity in the derivative of the scalar two-point integral: 

72 „ 1 . A'r 



l-Mfnl = -a-'^^g'^lM^l • (A4) 

The quark-squark loop with a color factor of 1 /2 contributes through two different fermion 
number flows due to the Majorana nature of gluinos. We take into account n/ = 5 light 
(s) quark flavors and a heavy top (s) quark. We do not include mixing in the top squark 
sector and take the top squark mass equal to the light squark masses. We set rrit = 175 
GeV. 

The gaugino couples only electroweakly to the quarks and squarks and thus does not give 
rise to strong self-energy corrections. All external particle self-energies have been renormal- 
ized on-shell and multiplied by a factor of 1/2 for proper wave function renormalization. 

The self-energy correction of the internal squark propagator depends on the off-shell 
squark four-momentum squared. Therefore, it factorizes only the corresponding t- or u- 
channel interference piece of the Born matrix element: 



^2 



— I 



^ — 772- ?77/- — 777,- 

Bo{t- m-g, 0)- 1 - Boiml m-g, 0)-^ / 

if: — m| H ^ -j. _ 

t + 2m? 
Bo{t; nig, 0)- ^ - Boinig; nig, 0) 



4CFgt[\M'M'*\ + |A^*A^"*|] 



t — 7n| ^ t — m| 



4Cf9;[\M'M'*\ + |A4*A4"*|]. (A5) 



The ultraviolet divergences cancel between the quark-gluino loop contribution and its su- 
persymmetric counterpart, the squark-gluon loop contribution. Since the gluon tadpole 
contribution is quadratic in the loop momentum it vanishes in dimensional regularization. 
The squark tadpole contribution vanishes after renormalization. The w-channel result can 
be obtained from the t-channel result given above through the exchange t ^ u. 

Like the self-energy correction of the squark propagator, the corrections to the quark- 
squark-gluino and quark-squark-gaugino vertices depend on the four-momentum squared 
of the squark and factorize only the t- or w-channel interference pieces of the Born matrix 
element. The quark-squark-gaugino vertex receives corrections through a gluon and a gluino 
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exchange between the initial state quark or antiquark and the squark that are proportional 
to the gauge and Yukawa coupling, respectively. For the t-channel we find 



|yW99x| = i [Bo{t; m-g, 0) - 5o(m|; m^-, 0) + Co(m|, 0, t; 0, m^, m^)(t + m? - - m|)] 



9 9X* 



t — X^^ ^ 



+ i 



f _|_ 777,^ ?72 - 

Bo{t; mg, 0)— \- - 5o(m|; m^-, 0)- ^ + Co(m|, 0, t; m,-, 0, 0) 



x(m| - mj-)] 4CFgl[\M'M'*\ + |A^*A^"*|]. (A6) 

Function Co{pl,pl, (^1+^2)^; "^1, "^2, "^3) stands for the scalar three-point integral defined in 
Appendix |B[ The ratio of quark-squark-gaugino couplings jX^^ ^ accounts for reversed 

fiavor flow in the vertex correction with respect to the underlying Born matrix element in the 
case of the exchange of an additional gluino. For neutralinos with real couplings it reduces 
to unity, whereas for charginos with real couplings it is given by a ratio of chargino mixing 
matrix elements. The infrared singularities 
^^2 



I A/f^'^^l 





'1 








- + 


-log 
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e 





1 2 

m~ — m% 
t — m| 



^CFg%M'M^*\ + (A7) 



arise from the gluon exchange correction. In dimensional regularization, the gauge bosons 
have (n — 2) degrees of freedom whereas their supersymmetric counterparts, the gauginos, 
have two. This difference leads to a mismatch between the quark-quark-gauge boson gauge 
couplings and the quark-squark-gaugino Yukawa couplings through finite next-to-leading or- 
der terms. The (super-)symmetry between the gauge and Yukawa couplings can be restored 
through a finite renormalization contribution 



-2 



l-^finLI =-C,CFgl\\M'M'^\\\M'M-% (A8) 

that can be found by comparing the quark-squark-gaugino vertex correction given above 
with the corresponding quark-quark-gauge boson vertex correction in exact supersymmetry. 
All M-channel results can again be obtained through the exchange t u. 

The quark-squark-gluino vertex correction can be obtained from the quark-squark- 
gaugino vertex correction if rriy. is replaced with rrig and the color factor Cp with Cp — Nc/2. 
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There are, however, two additional contributions with a color factor of Nc/'2^-, due to the 
non-Abelian gauge coupling of the gluino to the gluon, when a gluon is exchanged between 
the final state gluino and the initial state quark or antiquark and the squark. For the total 
quark-squark-gluino vertex correction in the t-channel we find 



\M 



<iqg\ 



+ i 



t — irt^ 



t — m? 



+ Co(m?,0,t; mg,0,0)(t-m2) 



Bo{t;0,m^ ^ 



2{t - m 



^ t — mi 



+ Co{m- 0,t■,0,mg,mg)- 



^N, 



t — 



(A9) 



It contains the following infrared singularities: 



\M 



=^2 



IR 



\MZ\\m^ ^ nig, Cf^Cf- Nc/2) 

12\ 1 



2^ + 7+2^^°^ 



m- 



log 



—t 

2 

mi 



(AlO) 



The finite renormalization contribution for the gluino vertex correction is 



\M 



mg I 

finite I 



C, -Nc - Cf gt[\M^M^*\ + \M^M^% 



(All) 



For the w-channel contribution, t and u have to be exchanged as before. 

Turning to the box diagrams that contribute to the associated production of gluinos 
and gauginos, we notice that they depend naturallly on the four-momentum squared of the 
exchanged squark. Therefore, we do not expect the full Born matrix element to factorize. 
In addition, the traces of the Dirac matrices project out terms that depend on the final 
state masses separately from those that do not, so that the squared t-channel and w-channel 
diagrams and the interference term can only be factorized individually. 

For the first box diagram in Fig. |^ we find 
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\M^^ = iBoit; m^, 0)it - ml)g^,iCF - Nc/2)4\WM^\/s 
+ iBoit; m,-, 0)(t - rnDgHCp - Nc/2)2\M'M'*\ 

{-2ml/{t-ml)/{t-ml)-2/{t-ml)) 
+ iBois; 0, 0){t - ml)/K{s, ml, ml)g^,{CF - Nc/2)2\M^M^\ 

{2s + A{t - ml) - 2m\ + 2m?) 
+ %B^{s- 0, 0)(t - m?)/i^(s, m% m^^gliC^ - NcI2)2\M}MF\ 

{%sm\l{t - m\)/{t - ml) + 2sm|/(i - m|) + 8sm?/(i - m|) 

-6sm?/(t - ml) + 2m|m?/(t - m|) - 2mJ/(t - m|) + 2m|m?/(t - m?) 

-2ml/{t-ml)) 
+ iBo{s; 0, 0)(t - ml)gl{CF - Nc/2)2\M'M'*\ 

{2ml/(t-ml)/(t-ml) + 2/(t-ml)) 
+ tBoiml; mg, 0){t - ml)/K{s, m|, m^^^^C^ - iVc/2)2pWVW^ 

{-2{t - ml)mys + 2{t - m?)m?/s - 2s - 2{t - ml) + 2m\ + 2m?) 
+ zBo(m?; m,-, 0)(t - m?)/ir(s, m|, m?)^2^C^ - iVc/2)2|A1 

(-4sm|/(i - m|)/(i - m?) - 2sm|/(t - m|) - 4sm?/(i - m|) 

+4sm?/(t - m?) - 2m|m?/(t - m|) + 2mJ/(t - m|)) 
+ iSo(m|; mg, 0)(i - m?)/X(s, m|, m?)g,'(Ci. - Arc/2)2|M*M«*| 

(2(t - m?)m|/s - 2(t - m?)/sm? - 2(t - m?) - 4m?) 
+ iSo(m|; m,-, 0)(i - m?)/X(s, m|, m?)^,'(Ci. - 7Vc/2)2|M*M**| 

(-4s/(t - m?)/(t - m?)m| - 4s/(t - m?)m? + 2s jit - ml)ml 

-2/{t - ml)mlml + 2/{t - m?)m|) 
- iBoiml; m,-, 0){t - m?)^2^C^ - Nc/2)A\M'M-*\/s 
+ iCo(m?, m|, s; 0, m,-, 0){t - ml)/K{s, m?, m?)^,'(CF - Nc/2)2\M'M^*\ 

(2m?s + 4m?(t - m?) - 2m?m? + 2m?m? + 2s(t - m?) + 2sm? 

-2(i - m?)m? - 2{t - ml)ml + 2m?m? - 2m|) 

39 



+ iCo(m?, m|, s; 0, rrig, 0){t - m1)/K{s, m|, ml)g^,{CF - Nc/2)2\M'M'*\ 
(8m?s/(t — m|) /{t — m^)m'^ + 2m?s/ (t — m^)m^ + Srn^s/ {t — m|)m? 
-6m?s/(i - m?)m? + 2m?/(i - m|)m|m? - 2m?/(i - m|)mj 
+2m?/(t - ml)mlml - 2ml/ {t ' ^IWg + 8s/(^ - ^)/(^ - ^D^I 
+6s/(i — m|)m|m| + 8s/(i — m^m- — 2s/{t — m^)rn'~.m^ — 8s/{t — m^)rn~ 
+2/ it - m|)m|m| - 2/{t - ml)m\ml - 2/{t - m^)m\m^- + 2/{t - m?)mjm?) 

+ iCo{ml m% s; 0, m,-, 0)(i - m?)^,2(CF - Nc/2)8\M^M^\ 

+ iCoiml ml, s; 0, m^-, 0)(t - m^gHCp - Nc/2)2\M'M'*\ 

i-mls/it - m\)/{t - ml) + 2m?/(i - m|)/(i - m?)m? + m?/(i - m|) 

-m?/(t - ml) + - m|)/(t - m?)m? - - m?) 

+2/(t - m|)/(t - m?)m| + 3/(t - m|)m? + m\/{t - ml) - 2/{t - ml)ml) 

- tCoiml, 0, t; mg, 0, 0){t - mDgliCp - Nc/2)4\WM^\ 
+ iCoiml, 0, t; mg, 0, 0){t - m?)^,'(CF - Nc/2)2\M'M'*\ 

{-l + ml/{t-ml)-ml/{t-ml)) 

- iCoiml, 0, t; mg, 0, 0){t - ml)gl{CF - Nc/2)4\WM^\ 
+ zCoiml, 0, t; mg, 0, 0){t - mDg^iCp - Nc/2)2\M'M'*\ 

(-l + ml/(t-ml)-ml/(t-ml)) 
+ tCo{0, 0, s; 0, 0, 0)(t - m?)^2(^^ _ iv^/2)2|>l*7W**| 

(mis/it - ml)/(t - ml) - s/(t - m|)/(i - m?)m? + s/(t - ml)) 
+ iDo{ml, ml, 0, 0; 0, mg, 0, 0)(t - m^^gHCp - Nc /2)A\JWM^\s 
+ iL'o(m|, ml, 0, 0; 0, m,-, 0, 0)(i - m?)^^^^^^ _ Nc/2)A\M'M'*\s 
+ zZ^o(m|, m|, 0, 0; 0, mg, 0, 0)(t - mfj^gHCp - Nc /2)2\M'M'*\ 

(-s/(t-ml)-s/(t-ml)) 
+ zDoiml, ml, 0, 0; 0, m,-, 0, 0)(t - m^^^^^C^ - iVc/2)2pWVW^ 

(./(i-m|)/(i-m?)). (A12) 
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Function K{s, m|, m?) = s^— 2s(m|+m?) + (m|— m?)^ is the triangle (Kallen) function of the 
partonic center-of-mass energy and the masses of the produced particles. The corresponding 
M-channel contribution is obtained by exchanging t and u. 
For the second box diagram in Fig. ^ we find 

\M^\^ = iBo{t; m-g, 0){u - m^^^^Cf - Nc/2)A\MKM^\/s 
+ iBo{t- m-g, 0)(m - m^^gl^Cp - Nc/2)2]WM^\ 

(2/ {t - w?^)/{u - m^^m\ + 2/{t - m?)/ {u - m^)m^g 

+2/{u - m\)l{u - ml)ml + 2/{u - m|)) 
+ iBo{s-m^,mg){u - mj^ / K{s,m\,ml)gl{CF - Nc/2)2\MKM^\ 

{-2s - 4(m - m?) + 2m\ - 2m?) 
+ iBo{s-m^,mg){u - mj;) / K{s,m%ml)gl{CF - Nc /2)2jM^^M^\ 

(—8s/ [u — m^) /[u — rngjuig — 2s/[u — m^)m^ — 8s/ [u — m^jrifig 

+6s/ (m — m?)m? — 2/{u — m|)m|m? + 2/(n — rn^^m^ — 2/{u — m'^)m'^m^ 

+2/{u-ml)ml) 
+ iBo{s- Trig, mg){u - m?)^^^^^ _ Nc/2)2\M^M^*\ 

{-2/{u - ml)/{u - ml)ml - 2/{u - m|)) 
+ tBoiml; mg, 0){u - ml)/K{s, m|, rnDgliCp - Nc/2)2\M'M^*\ 

{2{u - ml)mys - 2{u - ml)/sml + 2s + 2{u - m?) - 2m| - 2m?) 
+ z5o(m^; m,-, Q){u - ml)/K{s, m|, rnDgliCp - Nc/2)2\M^M^*\ 

(4s/ (m — m|) /{u — m^)ni^ + 2s/{u — m|)m| + 4s/('U — m^m^g 

— 4s/('U — nigjiTLg + 2/(-u — m^jm^nig — 2/[u — m^)m^) 
- iBoiml; rUg, 0){u - ml)g^,{CF - Nc/2)4\M'M'^*\/s 
+ tBo{ml; Trig, 0){u - m?)^^^^^ _ Nc/2)2\M^M^*\ 

(-2/(t-m|)/(M-m|)m|) 
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+ iSo(m|; m,-, 0)(t^ - m?)/X(s, m|, ml)g^,{CF - A^c/2)2|M*M"*| 
{-2{u - mfjmys + 2(m - m^g)/sml + 2(m - m?) + 4m?) 

+ iSo(m|; m^, 0)iu - ml)/K{s, m|, m?)^,2(C^ - 7Vc/2)2|At«M-| 
(4s/ (m — m|) /{u — m^)m^ + 4s/ (m — m'^)m^ — 2s/{u — m^)m^ 
+2/{u - ml)mlml - 2/{u - m?)m|) 

+ z5o(m|; rrig, 0){u - ml)gl{CF - Nc/2)2\M^M^*\ 
{-2/{t-ml)/{u-ml)ml) 

+ iCo(m?, m|, s; m^-, 0, mg-)(M - m?)/i^(s, m|, m?)5,^(Cj. - iVc/2)2pWVW^ 
(2m^s + Am^q{u — m|) — 2'm'^m'^ + 2m|m| — 2s{u — mf) — 2sm~ 
+2{u — m?)m| + 2(m — m?)m? — 2m|m? + 2m|) 

+ iCo(m?, m|, s; m,-, 0, m,-)(ii - m?)/X(s, m|, m?)^,2(Ci. - 7Vc/2)2|A^"A4"*| 
(8m?s/ {u — m|) /(m — m^)m'^ + 2m|s/ (m — m|)m| + Sm^s/ [u — m^^m^ 
—6m'^s/{u — m^)m^ + 2m?/(w — m'^)m^rn^ — 2m|/ (li — rn^)rn^ 
+2m?/ (m — m^)m^m^ — 2m?/(M — m^)m'^ — 8s/{u — m|) /{u — m'^)m^ 
—6s/ (li — rri^)rn^rn^ — 8s/ {u — m|)m| + 2s/ (i* — m^)rn^m^ + 8s/ (i* — m?)m| 
—2/ (m — m^)rn^rn^ + 2/{u — m|)m|m? + 2/(m — m?)m|m| — 2/{u — m^)rn^m^) 

+ iCo(m?, m|, s; m,, 0, mg)(i^ - m?)^2(Ci. - A^c/2)2|M*M«*| 
(1 - 2m?/s + 2(m - m?)/s - m|/s + 3m?/s) 

+ iCo(m?, m|, s; m,, 0, m^){u - rnDg^^iCp - Nc/2)2\M^M^*\ 

[—2m^s/ [u — m^)/{u — rrig) + Am^/[u — m^)/{u — mgjrrig + 2m^/ (m — m^) 
-2m~/ (li - m~) — s/{u — m|) /(x* - m^^)w?^ — s/{u — m|) + s/{u — m|) 
+3^ /{u — m|)/ (w — m?) — 4/(m — m'^)/{u — m'^)m'^ — 3/{u — m|)m? 
+3/(u - m|)m?) 

+ iCo(m2, 0, t; 0, m,-, mg){u - rnDgHCp - Nc/2)2\M'M^*\ 
(-1 + 2m?/s -{u- ml)/s - 2ml/ s) 
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+ iCoiml 0, t; 0, rUg, m~g){u - m?)^,'(CF - Arc/2)2|M«At«*| 

(2m?s/(M - m^)/{u - m?) + 2m?/(t - - m|)m| + 2m?/(M - m|) 

-s/ - m|)/ {u - m^g)m^g - s/{u- m|) - / {u - m|)/ {u - m^g) 

-2/ (t - m|)/ (m - m\)m\m'l - my {u - m|)) 
+ iCo{m% 0, 0, m,, mg)(ii - m?)^,'(Ci. - Nc/2)2\MM^\ 

( — 1 + 2m?/s — (m — w?g)/ s + m|/s — 3m?/s) 
+ iCo(m|, 0, 0, m,, mg){u - m^g)gl{CF - Nc/2)2\M^M^*\ 

(2m?s/ {u - m|) /(« - m?) + 2m?/ (t - m?) /(« - m^g)m^g + 2m?/ - m?) 

—s/{u — rn-) /{u — mf)ml — s/{u — ni-) — s^/ {u — m|) /{u — m|) 

-2/{t - mfj/ {u - mfim^g - my{u - m?)) 
+ iCo(0, 0, s; rUg, m~g, mg){u - m\^gl{Cp - Ncl2)2\M'M^*\ 
+ iCo(0, 0, s; mg, m^, mg-)(M - vr?^gl{CF - iVc/2)2|A<«M«*| 

(-2m|s/ (li - m|)/ (li - m?) + 5/(1* - m|)/ (i* - m?)m? + s^/ (li - m|)/ {u - m?)) 
+ iDo(m|, m?, 0, 0; m^, 0, m,-, m5)(M - m2)52^CF - iVc/2)2|A<*M"*| 

(— 2m^(ii — m^)/s + m^m^l s — hm~m^/ s — 2m^ + 2m^/s + 2[u — m^)/sm^ 

—m^/srn^ + 3/sm| + s + (m — m?) — m| + 3m?) 
+ iL>o(m|, m?, 0, 0; m^-, 0, m,, mg)(ii - ml)gl{CF - A^c/2)2|M«M"*| 

(— 2m?s/ (m — m?) /{u — m^)rn?g — m?s/ (m — m?) — 3m?s/ {u — m?) 

-4m?s^/ (li - m?)/ - m?) + 4m?/ {u - m'^)/{u - ml)m~ + my {u - m^^m^g 

—2>my{u — m?)m? + Am^s/{u — m^j/ {u — m?) — 2m|/(-u — m^)/{u — m^)m^g 

+2m|/ {u — m?) — 2s/{u — m|) /{u — m^)m'^ — s/{u — m^m^ + 2s/{u — m^)m^ 
/ {u-m'^)/{u- m^)m^g + /{u- m"^) + s^/ {u - m?) /{u- m"^) 

-2/{u - m\)/{u - m?)m| - m|/(ii - m|) + m|/(ii - m?)). (A13) 

The corresponding li-channel contribution is obtained by exchanging t and u. 
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For the third box diagram in Fig.^ we find 

JM^'^ = iBo{t; TUg, 0)(t - ml)g'^^Nc2\WW*\/s 
+ iBo{t- nig, 0)(t - m^^glNc2\M'M'*\ 

{-my{{t - m|)(t - ml)) - l/(t - m\)) 
+ iBoiu] m-g, 0)(t - m^^glNc2\M*M''*\ls 
+ i5o(«; m-g, 0)(t - m^^glNc2\M^M^*\ 

(m2/((t - m\){t - ml)) + m|/((t - m|)(M - m|)) + l/{t - m|) 

+m2/((t-m2)(M-m2))) 

- i5o(m?; mg, 0)(t - m?)^2jY^2|A^*A^«*|/s 
+ iSo(m?; m,-, 0)(t - m''^glNc2\WM^\ 

(-m|/((t-m|)(M-m|))) 

- iSo(m|; m^, 0)(t - m'^glNc2\M'M^*\/ s 

+ zSo(m|; m-g, 0)(t - m2)(72^c2|yW*^**|(-m2/((t - m|)(M - m^))) 

+ zCo(m?, 0, t- m-g, 0, 0)(t - m?)^2^c2|A<*A<"*| 
(-1 - (t - m|)/(2.) + m|/(2.) - m2/(2s)) 

+ zCo(m?, 0, t; m^-, 0, 0)(t - m?)^2^c2|A<*A<**| 
(-1/2 + m|/(2(t - m|)) - m|/(2(t - m^))) 

+ 2Co(m?, 0, u- 0, m,-, m^)(t - m^^glNc2\}^^M^\ 
(-1/2 + m|/s - (t - m^)/(2s) - m|/s) 

+ zCo(m?, 0, u; 0, m,-, m5)(t - m'-^glNc2\M'M'*\ 

(1/2 + m2s/(2(t - m|)(t - m^)) + m\my{{t - m|)(^i - m|)) 
+l/2m2/(t - m|) - sm2/(2(t - m\){t - m^)) + s/(2(t - m^)) 
-m|m2/((t - m|)(M - m|)) - m2/(2(t - m|))) 
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+ iCo(m|, 0, u] nig, 0, 0)(t - mf)glNc2\M'W*\ 

(1/2 + {t- ml) /{2s) - my {2s) + m|/(2s)) 
+ 2Co(m|, 0, m-g, 0, 0)(t - m2)^2^^2|^*7W**| 

(-1/2 - m5-./(2(t - m|)(t - m|)) - mj7(2(t - m^)) 

+sm|/(2(t - m|)(t - ml)) - s/{2{t - mj)) + m|/(2(t - m|))) 
+ ^Co(m|, 0, t; 0, m^, m,-)(t - m?)(72Arc;2|7W*7U«*| 

(-m|/s + (t-m2)/(2s)+m2/s) 
+ tCo{ml, 0, t; 0, m^, m,-)(t - m|)(7,'Arc2pWVW*^ 

(-1/2 + m?m?/((t - m|)(t - m?)) + m?/(2(t - m|)) - mj/((t - m|)(t - m?)) 

-ml/{2{t-ml))) 
+ iDo{m% 0, m?, 0; 0, m^, m^-, 0)(t - m''^glNc2\M'M'-*\ 

(-S- (t-m?) +m| -m?) 
+ iL>o(m|, 0, m% 0; 0, m^, m^-, 0)(t - m?)(72Art;2pWW^ 

(-S- (t-m?) +m| -m?) 
+ iDo{m% 0, m?, 0; 0, m^, m^, 0)(t - mf,^glNc2\WM^\ 

(_l/2 - (t - m|)/(2.) + m\/{2s) - m|/(25)) 
+ iL'o(m|, 0, m?, 0; 0, m^, m^, 0)(t - m?)2(72Ar^2|A^*A^**| 

(1 + s/{2{t - m|)) + s/{2{t -ml))- m|/(2(t - m|)) + m|/(2(t - m|))) 
+ iDo{ml, 0, m|, 0; 0, m-g, m^, 0){t - ml)^g^^Nc2\M'M'*\ 

{-s/{2{t - m|)(t - m|)) - l/(2(t - m|))). (A14) 

The corresponding w-channel contribution is obtained by exchanging t and u. 

The coefficients of the infrared singularities in the gluon exchange box diagrams can be 
simphfied considerably. The gluino-exchange box diagram is infrared finite. For the first 
box diagram in Fig. |^ in the t-channel we find 
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and the w-channel contribution can be obtained by exchanging t and u. The result is 



completely symmetric under the exchange rrig 



m 



The infrared singular pieces of the third t-channel box diagram in Fig. ^ are 



-a 



m?- — t 

2" 



1 f m~ —u 
f ^--log^ 



xA-fg',[\M'M'*\ 



|A^*A^«*|]. 



X 



mi 



(A16) 



For the w-channel result t and u have to be exchanged. 

Finally we sum the infrared singularities encountered in the virtual corrections and sep- 
arate them into Cp and Nc color classes: 



\M 



y I 

ml 



-a 



^ + ^ + ^log(^ 



ACFgi\M^\ 



— + - log '—^ 

e e \ smi 



(A17) 
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APPENDIX B: SCALAR INTEGRALS 

The tensor integrals that occur in virtual loop diagrams can be reduced to one-, two-, 
three-, and four-point integrals that are scalar functions of the loop momentum q The 
general scalar one-point integral is defined as 

The general scalar two-point integral is 

f cl"o 1 
5o(/;mi,m2) = / j^r-, 2^Tr—^, 2T' (B2) 



and its derivative is 



(27r)" [g2 -m2][(g + p)2 -m2]' 



d 



So(/;mi,m2) = ^So(g^; mi, ms' 



(B3) 



The general scalar three-point integral is 



^ f ^ 1 .g^. 

^ J (27r)"[g2-mf][(g + pi)2-mi][(g + pi+p2)'-mi]' 

and the general scalar four-point integral is 

X I . (B5) 

J (27r)" [g2 - ml] [{q + pi)^ - m|] [{q + pi+ P2Y - m^] [(g + P1+P2+ PzY - ml] 

The Pi,i = 1...4, are the four-momenta of the external particles, and the mi,i = 1...4, are 

the masses of the adjacent internal particles. 

The scalar two- and three-point integrals relevant for the associated production of gluinos 

and gauginos were calculated previously in a different physical context . We recalculated 

them and checked that the results in Ref. are correct. The divergent four-point integrals 



that contribute to the gluon exchange box diagrams in Fig. |^ were unknown, and they are 
presented here for the first time. We calculate the absorptive parts with Cutkosky cutting 
rules and the real parts with dispersion techniques. 
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The exchange of a massless gluon between the initial state quark and antiquark in the 
first box diagram of Fig. |^ leads to the following divergent four-point integral: 



Do{0, 0, m? m? s, t; 0, 0, 0, rriq) = iC^ 



s{t - m?) 



1 1 



In 



s(t — m? 



2\2 



(m? - m?)(m? - m^Q^ 



-2Li2 1 + 



2 2 

m| — m| 
t — m?~ 



- 2Li2 1 + 



2 2 

— 'm% 
t — m?~ 



Li2 1 + 



sm- 



TT 

-In 



In^ 



s 

2 



In^ 



In 



2 2 

m| — m| 



+ 2 In 
-In 



s 



In 



t — mi 



1 2 

m| — m| 



In 



2 2 

m| — m|, 

2 



(B6) 



This integral also contributes to the pair production of gauginos of unequal mass |I0|. In 
the limit of two final state particles of equal mass our result agrees with Ref. |^ . 



In the third box diagram of Fig. the exchange of a massless gluon between the the 
final state gluino and the initial state antiquark gives rise to a second divergent four-point 
integral: 



Do(0, 0, m|, t, u\ rriq, m-g, 0, 0) = iC^ 



-fin 



-u + m| 
rrigQ 



7r 



— + ln' 



+Li2 1 + 



u — m 



9 



2 2 

m| — m| 



-U + TTl; 
b 

rrigQ 

^^2(1+ 2/2 ZX 



3 

2 In 



-t + m? 

2 2 

m-: — m% 



-M -|- mi 

h 

rrigQ 



In 



-t + m? 



2 2 

y A, 



- 2Li2 1 + 



2 2 

— 

y A 
t — m? 



■ (B7) 



Again, our result agrees with Ref. for the case of two final state particles of equal mass. 

The supersymmetric counterpart of the first box diagram is the second box diagram 
in Fig. 1^. Since supersymmetry is broken and the gluino propagator and the two squark 
propagators are massive, this diagram does not have infrared divergences. The corresponding 
four-point integral can be expressed in terms of 16 dilogarithms [4^ . 
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APPENDIX C: SOFT GLUON EMISSION CONTRIBUTION 



In this Appendix, we collect the expressions for the finite pieces of the soft gluon emission 
contribution. The finite pieces remain after mass-factorization and cancellation of soft poles 
between soft gluon emission and virtual contributions, as described in the text. We begin 
with the Cp color class. 



dt2 du2 dt2 du2 V ^ 



Li2 
+ log 



[S + t2){s + U2) 
is + t2){s + U2) 



(CI) 



{■S + t2){s + U2 



sm 



log 



sm 



where 5 = A/m^ is the cut-off between hard and soft emission mentioned in the text. 
Its appearance in these terms is the explicit cut-off dependence that matches the implicit 
logarithmic behavior of the hard real emission contributions. 
The soft terms associated with the Nc color class are 



■2-5 



Ncas 

dt2 du2 dt2 du2 V 27r 



-2 + Lis 



U2t2 — s m\ 

{S -\- t2){s -\- U2) 



(C2) 



+ 



log 



{s + t2){s + U2) 



- 1 



log 



{s + t2){s + U2) 



Again logarithmic dependence on the hard/soft cut-off 5 is apparent. 
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APPENDIX D: HARD GLUON EMISSION CONTRIBUTION 



In this Appendix, we collect explicit expressions for the thirty-six matrix elements that 
contribute to the real emission of a gluon in the 2-to-3 partonic subprocess q q ^ g g X- 
The hard gluon emission cross section (in four dimensions) is 

dsidt2du2~ ds4dt2du2 IQn^ 36 52(54 + 771?) ^ ^ ^ ^ 

j=1..8 j=i..8 

where i and j label the diagrams in Fig. ^. The remainder of the factorization process is 
d^af _ Cf as asS{s + t2 + ui- S4) / fi'^{s4 + mf) \ ^^^^ 



ds4, dt2 du2 36 TT \^ si 

'sl-2Si{s + U2)+2{s + U2f\ f Xtt 



54(5 + ^2) / \{t-miy 

2 Xtu S mi 7772 Xu U2 [^4 U2 - Ui {s + U2)] 



+ 



{t - m|)[(A„ - s - t2){s + U2) + 334] [(A„ - s- t2){s + U2) + 554]^ 

sl-2Si{s + t2)+2{s + t2)^\ ( Xtt2[Sit2-h{s + t2)] 



■S4(s + t2) / \[{^t - S - U2){s -\- t2) + S : 



12 



2 Xtu S mi 7772 Xu U2 

"T/ o\r/A \/,,\, n~r 



iu - 777? ) [(At - S - 772)(S + ^2) + S S4] {u - 777? ) 2 



The elements Mf- are 



,>9 _ -lQCFT^Xtt2 (54 + ^772) 

(t - 777? )2 (5 + 7x2) 



Mf2 



8(C^-iVc/2)X, 

(t- 777?)(At-S- 772) 
X |s [-2 t2 (S + M2) + S4 (S + t2 + M2)] / + - ^2] / 

+ [At s (At + ti) + t2 (s {ml - 77i2) + At {t2 + U2- At) - 7x2^2)] / ( 77^ 

+ [At (S + 772 - t2) + t2 (At - 5 - 772)] / + [5 - t2] / (^) } 
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8NcXtt2{s + U2-ml^ 



S4 {t - m|) 



2,2 Al), 



(t-m|)2 J [ \u'UjA 



+ [2{2t2 + 2ml-s-U2)]i (—] - 21 ( — 

\U7aJ \U7A 

-16{CF-Nc/2)Xt.,s^m,m2 1 \ 
{t-ml){u-m1J \t'u'J' 

-16 CpXtu mi m2 \ 



S-to] I 



U' Uqa 



It 

277(34 + mf){Au — s — t2) ^2 7r{s4 + ml 



S4 [{Au - s -t2)(s + U2) + S S4] S4 (s + U2) 

+ [2*2]/ ^ 



(8C.X,){[A,(A, + «l/(4^) -[2A, + «jf (^) +A,/ 



MtA/ V S4 / V[(Aw - S - ^2)(s + 1*2) + 

2 At + ^ 1 



2 



[( A„ - s - ^2) (s + 1*2) + s S4] (5 + ^2) 
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\s4t-ml) 



Mis 



I [At{{s + t2){At + 777.2 - ^l) + "54^2) + ^2(554 - - " SU2 - t2U2)] I ( J 

L \t uyaJ 

+ [s + 12] i (^) + [At{ml + ml + S4) - t2is + U2- At)] I (^^^ 

-[t + ml + s4i{l)^, 

- [4t2{2At + ml-ml)]i (-^] - 2 [^2 + 2r77^] / ( — 

\t'U7Aj \U7, 



2 7r(At — S — 772)(S4 + 77if) 2 7r(s4 + 77if 

S4 [(At - s - 772) (s + ^2) + s S4] S4 (s + ^2) 

= [ iA,-s-U2)iA^-s-t2) ) - ' [f 

+ [A„ - t2] i ( ) + [A„ + At - s - 12 - 772] / ( — - — ) 1 , 

yu'UeAj \U6AU7aJ } 

MI7 = / ^ATc Xt, 7711 7772 \ r ^^^^ ^ ^2 _ ^ ^^^^ ^ _ _ ^^^^] I ( ^_ ) 
V S4 (77 - 77l|J J [ \t'U7Aj 



U7A 



+ [s + U2- 2At] — 2^(1) } , 



V (^i-"^l)(A«-S-^2 / I \f'U7Aj 



+2 [5 + ^2]/ + [-2At(s + i2) + 2A„7i2 + 2s7ii]/ 

V^*6A/ \UqaU7A 



7i7A 



+ [777^ + 777? - S + ^2 - 1*2] /(I) } , 



52 



SS4- St2 + SU2-t2U2]I 



Mit-ml){u-mlJj { \t 

iC/ra f 4:NcXtumim2 \ f / , ^. 2 
Ml,^^—-^^—-^j^[A^(s + U2) + s - 

+[s + t2- 2A„] / (—] + 2 /(I) 1 , 

^ _ 16 A^c s mi m2 (^4 + 2 m^) 
sl{t-ml){u-ml) 

^ . f -4A^cX,mim. \ r _ ^ ^ _s,-ml- 2,m\)] i (^) 

\s4[t-ml){u-mlJj { \u6aJ 

+ [s + t2]i( — ]+[u2-s]i{l)], 
\U6aJ } 

ML - ( (^^§^1^) { - -(^. - - ^ ^ - K^) } 

= — 2T7:- — <^ 2 [{t2 - A„ (s + U2) + sti]I [ 

V [t-ml){At - S-U2) J [ yWudAj 

+2[s + U2\i ( —] +[-2A^{s + U2) + 2sh + 2Att2\I ^ 



^ = ( I [t2 {U2 -A,) + A,s-s (.4 + + 3mD] / (—) 

Wit - ml) {u- ml) J { \u7aJ 



+ [s + U2]l{^]+[t2-s]I{l)}, 



UeAj \U7A 

-IGCpTi Xu U2 (s4 + mf) 
{u-ml)Hs + t2) 
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S{Cf-Nc/2)X^ 



X |[A^(s - U2) + A„(sMi + ^2(^2 + M2)) + sM2(m2 - m\) - tiuj] I 
+ [s- U2] i (^) + [A^{s + h) - u,{s + h)] i (—) 



8 Nc Xu U2 [s + t2- mj] 



A{Cf-Nc/2)X,,U2 



(u — m? Y 



-21 



t' 



^ !^4ui[s + t2]i (^-^^^ + 2[3ui + 2ml- 



(8 Cf Xu) I An[A. + «i] / f ^V) - [2 A„ + u,] I 



27r(s4 + mf)\ / A„(A„ + Mi)(s + M2) 



S4 J \[{Ai, - S - t2){s + U2) + S S4] 

2Au + ui ^1 



[(Ay - s - i2)(s + 1*2) + S + li2, 

Vs4(M-m|jy I \ u' J 

+ [Al{s + U2) + A^{s {ml - ml) 

+ U2 {ml - m\)) - t2U2)\ i i—;^ — ) 



+ [A„(s + t + u)+ U2{Au -s- ti)] i 



1 



U6A 



+ 4ii2[-2A„ - + m?] / ( ) - 2 [^2 + 2m^] / ( — ) I , 
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+ [s4(s + ^2) + At,(A„ - - 5m^ - 3^2) 

+ (A„ - s - i2)(-A„ + + 3m^ + s + ^2 + ^2)] / 



M6A 



+ 2[ii2 + m^ + m? -2A„] / f — ) +2/(1)1 . 

\UeAj ) 
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APPENDIX E: (ANTI-) QUARK EMISSION CONTRIBUTIONS 



In this Appendix we present the matrix elements for the contributions from emission 
of an additional light quark or anti-quark in the final state. These matrix elements may be 
obtained from those for gluon emission upon crossing one of the initial state partons and 
the final state gluon. However, we compute them ab initio in order to provide a check on 
the vahdity of the results. In all cases, the result of the crossing agrees with the exphcit 
computation. We limit our presentation to the quark emission contributions since the anti- 
quark emission contributions may be obtained from these expressions by replacements of 
t -(-^ u, t' ■(-^ u', and My everywhere. 

The quark emission contribution has two parts, consisting of the remainder of terms that 
are coUinear singular, and thus have the angular integrations done analytically, as well as 
a large class of terms which are coUinear finite and for which the angular integrations are 
done numerically. The matrix elements with coUinear singularities may be written in a form 
very similar to the gluon emission cross section, with 

ds^dt2du2~ dsidt2du2 16 tt^ 9Qs'^{si + ml) 
The remainder of the factorization is 



d^al asOisCp 



dS4 dt2 dU2 48 TTS^ 



1 + log 



1 / 2sl - 2si{s + U2) + (g + U2) 
2V {S + U2Y 



Xtt2 ^ 2 Xtu s mi 1712 



{t - m|)2 {t- m|)[(s + U2) (Au-s- ^2) + s S4] 

Xu U2 {u2Si-Ui{s + U2)) 



(E2) 



[(s 1*2) (A„ - s - ^2) + s 54]^ , 
After partial-fractionation, the coUinearly-divergent pieces of the hard matrix elements are 

I \^^6aJ \u'U6aJ s4{s + U2){Au - S - t2) + S S4\^ 

27r(s4 + mf) 



SaH-S + U2){Au - S -t2) + SS4] 



_ I ' 16CpXt„smim2 ^ f f f 1 ^ 27r(s4 + ml) 



{t-m\) j\ \u'ugaJ S4[{s + U2){Au - s - 12) + S S4] 
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2,2 2 
^ + — STUi 



V (saA + ^|r2)(^-^|) / 

—SS — 2A„S4 + SS4 — 2A„t2 + 27712^2 — 2m?t2 + 5^2 + 2S4t2 + 2t2^2 — SM2] I ( I 

= f^^^) |[4(«2 - AJ(-A„ - + mD(-A„ - + + ^2)] / ( ^) 

V ^2 / I \u'S3aJ 

+2[-6A„ + 4m? - 4m^ + 4^2] / (^) +4 / (^^) 
+ [2(-A„ + - 3m^)(-A„ + U2) 

+ (-A„ + m\-m\ + M2)(-2A„ + m? - m2 + s - S4 + ^2 + M2) 
+(-A„ + m? - m^)(-2A„ + m? - + s - S4 + ^2 + 3x^2)] I { —\ 
+2[-6A„ + 4m? - 6m2 + s - S4 + t2 + 4m2]/(1) + 6/(s3a) 
+2[- A„ + m? - 3m^] / (^^^ + 2/(k') | 

+ (^-^^^) {-A„[(A„ -ml + m^)(2A„ - m? + m^ - s + S4 - ^2 + ^2) 

+ (-A„ + m? - m2 - M2)(-2Au + m? - m2 + s - S4 + ^2 + M2)] ( ^2 J 
+ [(A„ - m? + m2)(2Au - m? + m2 - s + S4 - ^2 + 1*2) 
+2A„(4A„ - 3m? + 3m2 -5 + 54-^2 + ^2) 

+(-A„ + m? - m2 - ii2)(-2Au + m? - m2 + s - S4 + t2 + 1*2)] I ( -7- — ) 
+4 / (^) + 2A„[-3A„ + 2m? + s - S4 + ^2] / 

+2[6A. - - . + - «J I (^) - 6/(1) - 2A. / (^) + 2 / (£) } 
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+ f :^^] {-^u[{K -mj + m2)(2A„ - mj + ml - s + S4 - t2 + U2) 

+(-A„ + ml-ml- U2){-2Au + ml - ml + s - S4 + t2 + U2)] I ( — J 

+2[-6A„ + Sml-3ml + s-s, + t2-U2]i[^)+Ai(^^^ 

+2A„[-3A„ + 2ml + s-Si + ^2] / ( — ) + 2[6A„ -2ml- s + 84- h]H^) 

\U6aJ 

-6/Ka) - 2Aj (—] + 2i{u')], 

\U6aJ } 

^ = f 8{CF-Nc/2)XA ^ ^^2_ ^2^^^^ + _ ^2 + m^).(- A„ + U2) 

+2(-A„ + ml- ml){-A^ + U2)(-A^ + s + t2 + U2)+ U2(-A^ + U2)(-A^ + s + t2 + U2) 
+ {Au -ml + ml){s + M2)(-A„ + s + ta + ^2)] / 

+2[-3A„ + ml-ml + t2 + 2u2] I (^) + 21 (^^^ 

+ [2AuAu - 2Aumi + 2Aum2 + A^s - mis + mis - 2Ayt2 + mlt2 - mlt2 - 3A„ii2 
+mlu2 - mlu2 + SU2 + 2t2U2 + 2^21*2] I ( — — ) 

\S3Ay 

+[-4A„ + 2ml -2ml- s + '^h + 3x^2]^!) + 2/(s3a) + [^2 + U2] I (^-^^ | 

f -8{Cf - Nc/2)Xu\ rr / A , 2 2 \ , A 

+ I I {[ss4{-Au + mi - m2 - U2) + A„'U2(A„ -5-^2) 

\ SU2 / 

-A„s(A„ - ml + ml + U2) + 2A„(-A„ + s + t2)(A„ -ml + ml + U2) 
+(s + U2){-Au + s + i2)(A„ + ui)] i ( ) + 2[-3A„ + - + ^2 - ^2] / f ^) 

+2 i + [-4A„A„ - 2A„m2 + A„s + 2A„mi - mis + mis + 

-SS4 + 2Aut2 - rnlt2 + 777,2^2 + s^2 - 3A„M2 + 2sm2 + 2^2^21 1 
+[8A„ - 2ml + 2ml-s- 2t2 + Su2\I{l) - ^I{uqa) 
+ [-2A„ + s + t2]i (—] + 2i{u') ] , 

\U6A/ ) 



1 

u'SsA 
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/ 4iVcX^\ r ^^^^ -ml + ml)ssi + S4M2(A„ - ^2) + 2(A„ - mj + ml)s4{-Au + ^ 

+(-A„ + ml- ml){-Au + 1*2) (s + 1*2) 

+ (-A„ + - ml){s + M2)(-A„ + 5 + ^2 + ^2)] / ( J 

+2[s + U2 — S4] i ( ) + [2A„m2 — SAuS + Sm^s — Smls + 2A„S4 — 2m\si 
\ u' / 

+2m2S4 + 777.1^2 — ^2^2 — AuM2 + 2m^U2 — 4m2M2 — S4'U2 — t2U2 — U2U2\ I 



+ [-2ml + 3s - 2s4 + M2]/(l) + [-2m^ - ^2 - ^2] / ( — ) 
+ ^^^^^^^^ { [2s4A„(-A„ + mi - 7712 - ^^2) + A„S4U2 + ss4(A„ -ml + ml + U2) 
+A„(A„ + ui){s + U2) + {Au-s- i2)(A„ + ui){s + U2)] i ( ) 
+2[s + ti2 — S4] / ( — -T" ) + [2Aum2 + A^s + m^s — vri^s — ss — 2AuS4 + SS4 
+mlt2 - m\t2 - st2 + 3A„M2 - 2sm2 - 2t2U2] i 
+ [-2ml - s + 2s4 - 3w2]/(l) -[s + ^2] I 

-4NcXiamirn2 



u 



U6A 



{t-ml){s + U2) 

[-AuS + ss + 2A„S4 - SS4 + St2 - 254^2 - A„M2 + SU2 + ^2^2] I 

+[2A^ + s-t2]i(—] -2/(l)j, 

8CFXtt2 \ f 1^^^ 27r(s4 + ml) 



U'Uqa 



{t - ml){t - ml) ) \ Si 

A{Cf - Nc/2)Xtt2S^A \ f 2^ 2^ ^ una 
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+2[-3A„ + ml-ml + 2s - S4 + Sts + Sms] / ( ) 

\S3aJ 

+ [4A^ - 2A„mi + 2A„m2 - 4A„s + m^s - mis + + 2A„S4 - SS4 - 6A„t2 + 2mit2 
—2m2t2 + 3st2 — 2s4t2 + 2^2 — 8A^iM2 + 2miU2 + 4^2 — 2m2U2 + 5su2 

-2S4U2 + Qt2U2] i {—^ ) I 

+ ( - 2/(1) + 4 / (^) 

+ [4A^ - 2A„mi + 2A„m2 - 4A„s + mfs - rajs + + 2A„S4 - SS4 - 6A„t2 
+2mit2 - 2m2t2 + 3st2 - 254^2 + 2t2t2 + SU2] I ( ) \ , 

\U UqaJ J 

M<i ( 8{Cf-Nc/2)X,S,a \ I 

+ [-SS4 + 5^2 - 2S4t2 - S4M2 + t2M2] I | ' 

A>9 ( ^^{Cp-Nc/2)Xt^m,m2 \ a _l . ^ V A ^ ^ . ^ ^1 f M 

^26 = 7T 2"^ S (-^« + t2 + ^^2)(-A„ + S + t2 + ti2) / 



+ + [-A„ + t2 + U2] I ( — ) + I 



M^r. = 5— < 2S - 2S4 + t2 + U2\I { 

+ [ss4 - 2s4(-A„ + s + ta + M2) + (s + M2)(-A„ + s + ta + M2)] / ( ) 1 

\u'ssaJ J 



l(t-m|)(t-m|)A ' V^'^r; ^ ' W J ^' 

A squark width F is included to regularize a possible on-shell squark pole. As discussed in 
the text, r serves only to regulate the divergence of these interference terms, and it is taken 
to be very small compared to all physical masses and momenta. 

The finite pieces of the quark emission terms are evaluated directly from Eq. (p8[) with 
the angular integrals integrated numerically. They may be expressed as 
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Ml = Ml^ = Mf, = Mf, = Ml, = MI3 = MI5 = = Mfg = 0, (E5) 



TTTi — ^ 2T2A P *32 - s S32 - 2 S4 S32 + 2 ^2 S32 - 2 S32 - w6 S32 + t S5 - i 

— 3 t' 7712 +t'u2 — S4 ^2 — ^2 m'} 

/ — 4 Ci? 777 S3A \ f„ 2 r. ^ ^ / ^ II 

. 2 I 2 2 1 (2 «32 - S S32 - 2 S4 S32 + 2 t2 S32 " 2 S32 7i - Tib S32 + t S5 - i 

— 3 t' 772.2 + ^' '^2 — S4 ^2 — ^2 , 

/8(Q, - Nc/2)XtS^^ , 

I — ri \ 2T2\ r ) "^32 S4 + S32 « «6 + S S32 «7 + S32 ^ Ui 

+2 S32 7^7 + ^2 li' 7*7 - t' U2 U7} , 

7 — -— TT^ — ] 2T^\ r I * *32 S4 + S32U Ue-2 S4 5327*7 - t u-j 

V (S + ^2) (S^A + ) ^7A ^2 / 

-771^ t' U-! - S32 t' 7*7 + S5 t' 7*7 - S4 ^2 t*7 " ^32 1*6 '"7} 

■ 2TT\ 7 TTT I * *32 ^4 + S32 77 7*6 - 2 S4 S327*7 - t U-j 

V^e (S3A + "^l.r^) 777A (S + t2) / 

-77*1 ^' '"7 - S32 t' 7*7 + S5 t' 7*7 - S4 ^2 ^7 " ^32 t*6 t*7} 



77*2 



77*2 



5-^; r \ {—S S4 — ZSt — 2S4U — It U — U Uq + I 7*7) , 

V7*6 (t* - 77*1 J 777A (s + ^2) / 

/ -4 Ci. Xj„ 77*1 77*2 S3A \ r , . r, I ' 1 O ' + ' 

, ^2r2W„ ^2 W {-SS32 + SS4 + 2SU - 2 S32 « + 2 S4 77 - ^2 « 

+2 7*'^ + 7*2 + 7*' 7*6 - t' 7*7! , 

/ -4 Xc ?7*i 77*2S3A \ 1 o / / , ^/ 1 

U.(^l. + m|.n(.-0(. + ..) j - 2..« - « «a + t «,} . 

/ 4Cj;'X„S4 7*2 \ fr. 2 r. 2 / 1 

{2-2 + 2mi-.32 + .4-.5-^.+^.2-M> 

-8 {Cf - Nc/2) Xtu S4A SsA mi 77*2 \ ^ _ , _ r, 2\ 
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^ {^m\ s + m\ s — s s^2 — s + t2u' + t' U2 



-2 ^2 1*2 + 2 li' 1^2 - 2 + U2 Ue + t2U7 + 2 U2 U'j^ , 

—2 u' U2 + u' UQ + t' ui + 2 u' ur} , 

f 2NcXtuS4AS3Amim2 \ r o 2 / « 2 / 

V - ^i) («1a + ^|r2) + m|r2) lie; 

+2 s li' — 2 S4 li' + 2 S5 li' + t2 — ^' U2 — u' UQ + t' ut + 2 u' u^} , 

f -2NcXuS4AU2 \ f r, 2 R 2 . 02,0 

7 2 \2 / 2 — \ 2TT{ — I -2 7712 S4 - 6 mi S4 + s S4 - 2 S4 + 2 S4 S5 

—ml t' — 3 ml t' — 2 s^t' + S5 1' + S4 ^2 — 532 uq — u' ug + 2 S4 U7 + t' tiy} , 

/ 4:CFXtSs2U7 \r 2,0 2, ' ^ \ 

„2 /,2 , ^,2^2^ {2m2 + 2mi + S32-S4-S5-li -1^2 + 1^7}, 



/ 4 (Cf - Nc/2) Xtu S4A S3A mi ms \ r , o ,/ / , , / , o '2 , +/ 



+2 Xi' U2 — u' Uq — t' — 2 u' U-j} , 



f -A{CF-Nc/2)Xts,2 \ r 2^0 2 ^ ' ^ 

+t' U-j + t2U-j + 2u' U-j + 2U2U-J — 2UqUt — 2 j , 

/ 2A^c-'^tS32 \ r 2 2 2 / 

-r^ 2~n5\ |S m2 + 3 S m^ - S S4 - S S5 + 2 m2 + M Me + M2 «6 

+6 u-j m\ + 2 sut — 2 S4U7 — 2 S5U7 — t' U7 — t2U7 + 2 u^} , 

f -2NcXtuS3Amim2 \ . , ' 2 , « / 2 o ' 

■ 2-F?\ ? 2^ |SS32 - SS4 + 2m m2 + 67* mi - 2S4li 

—2 S5 7i' — 2 t' u' — t2u' — t' 7*2 + 7*' 7*6 + + 2 7*' 7*7} , 
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Ml, 



Ml, 



Ml, 



M^s 



—8 Cp Xu S4 12 

16 Cp Xtu S4A S3A mi 1712 t' 
4 Nc Xtu S4A S3A mi 1712 

s {sl^ + m|r2) U6 (sIa + ^|r^) 



^ {— s S4 + 2 s t' + 2 s u + u uq + t' U7 + 2 u U7} , 

..q f 4iVcX„S4A \ r 2 2 , 

Mg^g = ^ 2"^ -m — \ 2~F?V ] {-sSiTn^- ssim^ + ss^s^- SS4U2- S4 1*2 Uq 



+u' U2 Uq — S4 ^2 M7 — 2 S4 M2 M7 — t' U2 M7} , 



Xt .S;-!2 Uq 



2 ]\l X \ 

, 2 ^' '2 T^2A 1-4 sm? + 4 sii6 + 4 5-^7 + 4x^6 1*7 + 4ii?}, 

A^Q ( 2NcXtuS3Amim2 \ . , , , o ' \ 

„ f -4:NcXtS32 \ f 2 2-1 



-16 Nc Xtu S3A mi 777.2 u' (2 + Me) 
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TABLES 



mi/2 (GeV) 



rug (GeV) 



(GeV) 



Tevatron (pb) 



LHC (pb) 



100 
100 
100 
100 
100 
100 
200 
200 
200 
200 
200 
200 
300 
300 
300 
300 
300 
300 
400 
400 
400 
400 
400 
400 



284.9 

284.9 

284.9 

284.9 

284.9 

284.9 

533.7 

533.7 

533.7 

533.7 

533.7 

533.7 

775.5 

775.5 

775.5 

775.5 

775.5 

775.5 

1012.4 

1012.4 

1012.4 

1012.4 

1012.4 

1012.4 



31.1 

62.6 
-211.3 
240.9 

56.8 
240.2 

77.6 
146.6 
-358.2 
380.6 
145.2 
379.5 
120.3 
232.3 
-510.9 
528.2 
231.7 
527.2 
162.2 
317.1 
-665.3 
679.2 
316.77 
678.4 



0.165 

0.233 
0.00588 
0.00424 

0.329 

0.0100 
0.00283 
0.00231 

1.22- 10-5 
9.73-10-6 

0.00238 
1.680-10-5 

5.51- 10-5 
2.07-10-5 
1.46-10-^ 
7.71-10-9 
1.58-10-5 
1.06-10-^ 
6.58-10-^ 
6.85-10-^ 

5.23- 10-13 
9.85-10-1^ 

3.52- 10-^ 
1.16-10-13 



2.82 
4.50 
0.176 
0.364 
9.61 
0.928 
0.212 
0.354 

0.00452 
0.0151 
0.720 
0.0379 
0.0400 
0.0605 

0.000385 

0.00144 
0.128 

0.00368 
0.0110 
0.0151 
5.53-10-5 

0.000207 
0.0327 

0.000550 
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TABLE I. Total cross sections in pb at the Tevatron and at the LHC for the parameters of 
the SUGRA model. Shown in column 1 are four chosen values of the common fermion mass. In 
column 2 are the derived values of the gluino masses, and in column 3 are the masses of the four 



neutralinos and the two charginos. See also Figs. 11 and 12. 
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FIGURES 




FIG. 1. Leading order Feynman diagrams for qq ^ gx- 




FIG. 2. Generic virtual one-loop diagrams for the associated production of a gluino and a 
gaugino. The crossed regions denote self-energy and vertex corrections that are present only one 
at a time at next-to-leading order. 




FIG. 3. Self-energy diagrams for external quarks (g) and gluinos (g) and internal squarks (g) 
including the full supersymmetric QCD particle spectrum. Contributions from the tadpole graphs 
vanish. 
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FIG. 4. Vertex corrections for the quark-squark-gluino (qqg) and quark-squark-gaugino {qqx) 
vertices including the full supersymmetric QCD particle spectrum. 
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FIG. 5. Box diagrams for associated production of a gluino and a gaugino. There are u-channel 
box diagrams in addition to the diagrams shown. 
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FIG. 6. Feynman diagrams for the gluon emission contribution, qq ^ g g X- 
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FIG. 8. Dependence of the total hadronic cross section on the cutoff parameter S for the gxl 
channel in the SUGRA model at the Tevatron, with rUg = 1012 GeV and m^o = 679 GeV. 
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FIG. 9. Scaling functions for the SUGRA channel, with rrig = 410 GeV and m^o = 104 GeV. 
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FIG. 10. Scaling functions for the gxf channel in the SUGRA model with rrig = 410 GeV and 
101 GeV. 
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FIG. 11. Predicted total cross sections at the Tevatron for all six gx channels in the SUGRA 
model as functions of the mass of the gluino. See also Table I. 
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FIG. 12. Predicted total cross sections at LHC energies for all six gx channels in the SUGRA 
model as functions of the mass of the gluino. See also Table I. 
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FIG. 13. Calculated enhancements of the cross sections from NLO contributions at the Tevatron 
for all six SUGRA gx channels as functions of the mass of the gluino. 
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FIG. 14. Calculated enhancements of the cross sections from NLO contributions at the LHC 
for all six SUGRA gx channels as functions of the mass of the gluino. 
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FIG. 16. Dependence of the predicted NLO and LO total cross sections at the Tevatron on the 
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FIG. 17. Dependence of the predicted NLO and LO total cross sections at the LHC on the 
renormahzation and factorization scale. We select gx2 production in the SUGRA model, with 



rrig = 410 GeV and m^o = 104 GeV. 
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FIG. 19. Predicted total cross sections at LHC energies for all six gx channels for a gluino with 
mass 30 GeV as functions of mi/2- 
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FIG. 20. Calculated enhancements of the cross sections from NLO contributions at the Tevatron 
for all six gx channels with a light gluino of mass 30 GeV as functions of mi/2- 
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FIG. 21. Calculated enhancements of the cross sections from NLO contributions at the LHC 
for all six gx channels with a light gluino of mass 30 GeV as functions of mi/2- 
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FIG. 22. Differential cross section in rapidity da/dy for tlie production of Xi witfi mass 101 
GeV in association witfi a ^ of mass 410 GeV at the Tevatron. 
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FIG. 23. Differential cross section in transverse momentum da/dpT for the production of Xi 
with mass 101 GeV in association with a 5 of mass 410 GeV at the Tevatron. 
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FIG. 24. Differential cross section in rapidity da/dy for tlie production of xf with mass 101 
GeV in association witfi a ^ of mass 410 at the LHC. 
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FIG. 25. Differential cross section in transverse momentum da/dpT for the production of x-^ 
with mass 101 GeV in association with a 5 of mass 410 GeV at the LHC. 
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FIG. 26. Differential cross section in transverse momentum da/dpr for the production of X4 
with mass 309 GeV in association with a ^ of mass 410 GeV at the LHC. The sohd hne is the full 
NLO cross section. The dashed curve shows the LO Born cross section. The dotted line shows the 
contribution of the qg initial state subprocess. The dot-dashed curve is the difference between the 
full NLO result and the sum of the Born and qg contributions. 
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